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SOME EXAMPLES OF GLOBAL INSTABILITY 
OF THE COMPETITIVE EQUILIBRIUM* 


By HERBERT SCARF’ 


IN THIS PAPER we shall consider the problem of stability of the 
competitive equilibrium. The market demand functions are sums of 
individual demand functions obtained directly by utility maximization. 
The rate of change of the price of each commodity is assumed to be 
proportional to the excess market demand for that commodity. A 
number of examples are given for which the motion of the prices is 
globally unstable in the sense that starting from any set of prices 
other than equilibrium, the prices oscillate without tending towards 
equilibrium. 


1. INTRODUCTION 


The problem of stability of the competitive equilibrium is described 
in [3], and we shall content ourselves with a review’. Several individ- 
uals with utility functions U,(x,,-++-,2,) for the same commodities 
are engaged in the trade of these commodities. Each individual begins 
the trading with an initial endowment of goods, say J}, Ii, --+, I‘, for 
the ith individual. An initial vector of prices p,, +++, p, is announced, 
and each consumer then determines his demand for all of the com- 
modities by the usual procedure of maximizing his utility function, 
subject to the constraint that his expenditure shall not exceed the 
value of his initial endowment of goods at the stated price vector. 

For each commodity the sum of the individuals’ demand functions 

* Manuscript received January 21, 1960. 

1 Research undertaken by Stanford University under Contract Nonr-225(28), NR 047-019 
with the Office of Naval Research and completed at the Cowles Founcation for Research 


in Economics under Contract Nonr-358(01), NR 047-066 with the Office of Naval Research. 
I would like to thank K. J. Arrow, L. Hurwicz, and H. Uzawa for a number of stim- 


ulating conversations. 
2 The reader may wish to consult the excellent bibliography given in [3] and also in [1]. 
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minus the sum of the initial endowments of that commodity is called 
the market excess demand for the commodity; in this paper it will be 
denoted by f,(p,, +++, P,), Where the subscript refers to the commodity 
in question. (These functions are homogeneous of degree zero and 
satisfy the Walras Law }>p,f,=0.) Excess demand functions for 
individuals will always be denoted by 72,. 

An equilibrium price vector is, of course, a vector of prices for 
which all of the excess demand functions vanish, and recent work in 
this area [2, 8] has shown that under suitable regularity conditions 
equilibrium prices will always exist. (Because of the homogeneity of 
the demand functions, any positive multiple of an equilibrium is again 
an equilibrium, and in this sense we should speak of an equilibrium 
ray. There may be several such rays.) The stability problem, on the 
other hand, is less concerned with the existence of equilibrium and 
more with the question of what happens to the prices if initially they 
are different from the equilibrium. 

There is nothing in the model described so far which enables us to 
compute the motion of prices if we are not at an equilibrium point. 
For this we need some specific assumptions on the price adjustment 
process, that is, the procedure by which prices may be expected to 
change if we are away from an equilibrium price. The intuitive no- 
tion that an excess of supply over demand should result in a decrease 
in price, and an excess of demand should result in an increase in price, 
has been formalized mathematically in [7] and [3] by the statement 
that 


dp, _ eee 
“a = Al f(y. ’ Pn)) ’ 


where H, is a sufficiently regular sign-preserving function of its argu- 
ment. (These equations are to hold if all p,>0.) In order to be 
concrete, in this paper we shall generally take the functions H to be 
equal to their respective arguments, so that for each commodity the 
rate of change of price is equal to the excess demand. 

Now let us turn our attention to the problem of stability with this 
type of adjustment process. Early work in this area [5, 6] tends to 
emphasize what might be called ‘‘local’’ stability; the initial prices are 
assumed to be close to some equilibrium point, and an analysis is made 
of whether there is a tendency to converge to the equilibrium point, 
depart from the equilibrium point, or perhaps even a tendency to more 
complex types of behavior. The ‘“‘local’’ analysis proceeds by means 
of the linear terms of the Taylor series expansion about the equilibrium 
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point, thus converting the problem to a linear differential system with 
constant coefficients. It is possible in this type of analysis to obtain 
many examples of completely unstable equilibria (though it should be 
mentioned that there is an uncomfortable tendency for examples to be 
produced without any consideration of their origins as market demand 
functions derived by the summation of individual demand functions). 
The local analysis, however, is somewhat unsatisfactory, in that it is 
quite possible for other equilibrium points to exist, and the system 
cannot be said to be unstable without examining whether the prices 
tend to another equilibrium point. 

This consideration leads naturally to the problem of stability in the 
‘“‘global’’ sense [3], which is concerned with the solution of the differ- 
ential equations based on the price adjustment mechanism for an 
arbitrary initial set of prices. If the solution of the differential equa- 
tions approaches some equilibrium point as time becomes infinite, then 
we have global stability. Clearly it is quite possible for there to be 
several equilibrium points, none of which is completely stable from the 
local point of view (that is, attracts all neighboring points) and where 
the system, in its entirety, is globally stable. 

Thé problem of global stability, as it has been described above, has 
been discussed by several 2uthors. Though alternate methods of proof 
have been found, the results do not go essentially beyond those ob- 
tained in [1]. There are, first of all, several special cases which can 
be mentioned: global stability is known if there is a single consumer, 
or if there are only two goods, or if it happens that at some equi- 
librium point there is no trade. There is also an unpublished result 
of the author’s showing that global stability is known if all of the 
consumers have the same homogeneous utility function, but with 
different initial distributions of stock. Aside from these very special 
cases, the most important result is that global stability will occur if 
all of the goods are gross substitutes; mathematically, this means 
of,/ap,>0 if i# 7. Apart from this, very little else has been found; 
no assumptions substantially different from gross substitutability have 
been shown to imply stability, and up to the present no examples of 
instability have been produced. 

This paper presents a series of examples, all derived from utility 
maximization, which are globally unstable. The examples given here 
involve three consumers and three goods, but the techniques may be 
extended to either more consumers or more goods. The examples all 
involve a rather simple relationship between the utility functions of 
the three consumers, namely 
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(1) U2, La, £3) = U,(%,, Ls, 2) , 

U,( dy, La, Z) = UX, Xs, %,) , 
and a corresponding cyclic permutation of the initial endowments. 
Specifically, we will assume that if the initial endowments of the first 
consumer are (J,J, K) for the first, second, and third goods respec- 
tively, then those of the second consumer are (K,J,J) and those of 
the third consumer (J, K, I). This relationship is introduced to simplify 
the calculations, and, as the method will show, it may be relaxed. 

It is an easy consequence of utility maximization that for the price 
adjustment process we have selected (rate of change of price equal to 
excess demand) the motion of the prices will always be constrained to 
the sphere p? + p; + p} = constant. In all of our examples there will 
be an equilibrium point at (p, = p, = p,), which, except for the first 
example, will be locally completely unstable (there is a small region 
around this equilibrium point such that if the initial price is in the 
region, the prices will eventually leave this region and stay out). In 
the first example the motion of the prices will be in concentric curves 
about the equilibrium point, whereas the more complex examples will 
give rise to limit cycles. 

As we shall see from the examples of Section 3, instability does not 
depend on a delicate assignment of values of initial stocks or parame- 
ters in the utility functions. Though it is difficult to characterize 
precisely those markets which are unstable, it seems clear that in- 
stability is a relatively common phenomenon. 

What are the implications of these examples? It seems to me that 
there are several possible interpretations that might be made. 

1. One possible interpretation is that the model is substantially re- 
alistic and that instabilities of the type described in this paper could 
possibly occur. An even more presumptuous interpretation along these 
lines is that instability is responsible for some aspects of the business 
cycle, though for this sort of interpretation it would seem advisable 
to produce examples of instability with a model of a complete economy 
rather than a pure trade model alone, and this may be more difficult 
to do. 

2. Another possible interpretation is that the price adjustment pro- 
cess postulated above is not correct. This view can, of course, be heid 
without any reference to the question of stability. An argument for 
this position is that in one sense or another we are considering a 
dynamic process, and yet nowhere do the simplest dynamic considera- 
tions such as saving, interest, etc. appear in the model. 

8. As a final interpretation it might be argued that the types and 
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diversities of complementarities exhibited in this paper do not appear 
in reality, and that only relatively similar utility functions should be 
postulated, and also that some restrictions should be placed on the 
distribution of initial holdings. This view may be substantiated by 
the known fact that if all of the individuals are identical in both 
their utility functions and initial holdings, then global stability obtains. 


2. A VERY SIMPLE EXAMPLE OF INSTABILITY 


In this section we shall describe a very simple example which leads 
to instability. It will be seen that this example is quite an extreme 
one. In the next section we shall describe a number of additional 
examples which, while somewhat more complex, do not have the dis- 
agreeable features of the present one. 

Let the utility function of the first consumer be U,(x,, 2,, x;) = 
min (x,, 2.) and his initial endowments (1,0,0). This consumer has no 
desire for the third good and his indifference curves for the first two 
goods are of the form: 


Xe 














Xi 


FIGURE 1 


For any income M the same quantity will be demanded of goods one 
and two and therefore the demand functions are 





M 
Y:(Pi, Par Ps, M) = — , 
Dp, + Ds 
M 
YAD:, Po» Ps», M) = ———_-, 
DP, + Ds 


YAP, Po» Ps» M) = 0. 
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Now the income of the first consumer is derived from his endow- 
ment of a single unit of good one, so that M = pp, and therefore the 
excess demand functions of the first consumer (supply is being sub- 
tracted off now) are given by 








1,= —P , 
DP, + Ps 

t= P : 
Pi + Ds 

e=0. 


The excess demand functions of the second consumer are obtained 
by a cyclic permutation of all subscripts, 1— 2,2-—-3,3-—1, and once 
more for the third consumer. If we add these together, the market 
excess demand functions are given by 








f= —P; + Ps : 
‘D+ PtP 

h= a P; : 
’ DP; + Ds Py + Pp, 

f= —D, + Ds 





PtP B+’ 


and of course the price adjustment process leads to the differential 
equations 

d 
(2) SPs = f(D, Py» Ps) « 

dt 

It is a trivial matter to verify that p, = p, = p, is the only equilibrium 

point. The fact that }> pj = constant follows from }> pdp,/dt = 0, 
which is the Walras Law. In order to show that the solutions are 
unstable we shall demonstrate that p,p,.p, = constant for any solution 
of (2). This would follow if we could show that 


(3) S:DPs + f:PiPs + foD:Ps 
equals zero, but (8) is equal to 





PAPi — Pi) , PAPi— Pi) . Pi(Pi — Pi) 


P, + PD, Ps + Dy P, + Ds 
- DAP, sd P.) + p{p, — P,) + PP, = Ps) =0. 


That this implies instability is clear. Let the initial prices be chosen 
such that p? + p} + p; = 3, so that the intersection of the equilibrium 
ray and this sphere is (1,1,1). The value of p,p.p, at equilibrium is 
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one, and therefore if the initial price gives a value different from one 
to p,p,p, we never reach equilibrium. It should be remarked that the 
maximum of p,p,p, subject to the constraint )> pi = 8 is actually one, 
so that if the initial position is anything other than (1,1, 1), the path 
is completely unstable. 

Occasionally the following price adjustment process is discussed. 
One of the goods is singled out and its price kept constant; the re- 
maining goods are meant to vary according to the differential equa- 
tions given above. If, in our case, we put p,=1, then we are led 
to the system 





Oe oc cong 2d 
dt ptp, 1+p,’ 
dp, _ -—1 Pp; 








+ > 
dt P+1l BM+D, 
Routine calculations show that for this system 


p,p,e-2 i+? = const., 
and this is again sufficient to show instability. 


3. A CLASS OF EXAMPLES 


The example of the previous section has a number of special prop- 
erties. All of the Slutsky terms [dy,/0p, + y,(@y,/9M)] are zero, the 
indifference surfaces are not strictly convex, and certainly not differ- 
entiable, and, finally, the initial holdings are of a rather extreme type. 
It might be thought that one, or several, of these properties is re- 
sponsible for the instability of the previous example, and that stability 
would return if these properties were removed. The examples of the 
present section show, however, that it is quite easy to obtain instabil- 
ity with none of the objectionable properties mentioned above. 

We shall make an attempt to keep the reasoning relatively general 
in this section. Some conditions will be described which imply in- 
stability, and we shall demonstrate by specific examples that these 
conditions may be satisfied. 

As was mentioned in the introduction [see equation (1)], the utility 
functions of the three consumers will be obtained by a cyclic permuta- 
tion of the goods and the initial endowments. This means that if the 
excess demand functions of the first consumer are represented by 


L(Pi5 Poy Ps)y LkPry Pay Ds), Ta Pry Poy Ps), 


then the excess demand functions of the second consumer for the 
first, second, and third goods respectively will be given by 
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L( Po» Por Pr)> (Por Psr Pi)» LlPor Ds» Ps) » 
and those of the third consumer by 


LPs» Pry Pr), Cs Psy Pry Pr), (Ps Pry Ps) - 


(Continuity and differentiability properties of these functions will be 
assumed whenever necessary.) 

The market excess demand functions are obtained by the summation 
of individual demands and are therefore given by: 


F(Prs Pay Ps) = L(Dry Pay Ds) + LslPoy Pay D1) + LPs» Pry 2) » 
(4) SA Pr» Poy Ps) = LAPry Poy Ps) + L(Po» Psy P:) + Ls(Ps» Pry Pa) » 
SAPry Pos Ps) = LAD» Poy Ds) + (Pay Psy Pi) + L(Ds» Pr» Da) - 


The differential equations of stability are, of course, given by dp,/dt = 
S&Pv» Po» Ps), When all prices are non-negative. 

As we mentioned previously, the Walras Law (}> p,7, =0 for each 
consumer, and its consequence, )> p,f, = 0) implies that d(}> pi)/dt = 0, 
se that the motion of the prices will always be on that sphere with 
center at the origin which passes through the initial price vector. Let 
us assume that the initial prices have been selected such that }> pi=3. 
When we speak of equilibrium prices we shall mean the intersection 
of the equilibrium ray and this sphere. 


LEMMA 1. The price vector (1,1,1) is an equilibrium price. 


This is immediately obvious by an application of the Walras Law to 
equations (4), when all of the prices are set equal to one. 

It is not at all correct, for a general selection of the individual 
demand functions 2,, x,, #,, that the market demand functions f,, f,, f; 
(4) must have a unique equilibrium point at (1,1,1). There are a 
number of very simple conditions, however, which imply that this 
equilibrium is unique. We shall give only one such set of conditions. 


LEMMA 2. Let 
A = 6f,/Op, < 0, and B = 6f,/ap, > 0, 


everywhere in the orthant p, >0,p, >0,p, >0. Then the equilibrium 
point (1,1, 1) is unique (aside from positive multiples). 


In order to demonstrate this lemma, let us first make the observa- 
tion that 


FAP: Pay Ps) = F.(Par Ds» Ps) » 
FAPr» Pas Ds) = FAPa» Ps» Ps) 
FAP1 Par Ps) = FA Pr» Pr» P:) - 
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This implies that if (a@,8,¥) is an equilibrium point, then so is 
(8,7, a), and also (y,a@,). Therefore, if we have an equilibrium point 
different from (1, 1,1), then we may find an equilibrium point (a, 8, 1), 
perhaps by permutation, with either 


(5) a>B2Y7 
(at least one of the inqualities is strict), or else 
(6) asR<sY 


(again the same remark). 

Let us assume that the former holds. An argument similar to the 
one we are about to give works if instead of (5) we have (6). 

We shall show that f,(a, 8,7) must be different from zero. Now 


FAQ, B, Y) = v{a, B, Y) T 2,(8, 7, a) + x(Y, a, 8) . 


The assumptions of this lemma imply that as a function of a, f, is 
strictly increasing, and as a function of vy, f, is strictly decreasing. 
Since a > 8 >y7 (with inequality somewhere), we have 


FAQ, B, ) > 2(8, B, 8) ~ x,(8, B, 8) + x,(8, B, 8) =0 ’ 


by the Walras Law. This completes the proof of the lemma. 

Subsequently, in this paper, we shall have occasion to consider 
specific examples of demand functions for which the conditions of 
Lemma 2 do not hold for all prices in the orthant p, >0, p,>0, p, >0, 
but rather for certain subsets. We shall then have recourse to the 
following lemma, which involves a type of subset convenient for us. 
The lemma is demonstrated in the same manner as Lemma 2. It 
should be noted that the operations involved in the proof of Lemma 2 
do not take us out of the set described in Lemma 3. 


LEMMA 38. Let the conditions of Lemma 2 hold for the set of prices 
(p,/p; > &, P./P, > &, p,/P, > &), where ¢ is a small positive number. 
Then there are no equilibrium points in this set other than (1,1, 1) 
and its multiples. 


We shall now give a condition on the demand functions which 
implies that the equilibrium point (1,1, 1) is unstable. 

LemMa 4. If C= 4f,/@p, is positive, at the point (1,1, 1), then this 
equilibrium point is locally unstable in the following sense: there is 
a region (on the sphere >> pi = 3) about the point (1,1,1) such that 
any point in this region (other than the equilibrium point) will 
move away from the point (1, 1,1). 








166 HERBERT SCARF 


The Jacobian of the functions f,, f,, f,, at the point (1, 1, 1) is given by 
(° A 

BCA 

lupe 

all evaluated at that particular point. Moreover, the functions /,, f,, f; 


are homogeneous of degree zero, and therefore the Euler relationship. 
at the point (1,1, 1), tells us that 


(7) A+B+C=0. 


The Jacobian, therefore, has one characteristic root equal to zero. 
The other characteristic roots may be shown to have positive real 
parts, because of the assumption C>0. It is then possible to apply 
standard theorems of differential equations (see [4], chap. 13), in order 
to deduce local instability. However, the various properties of demand 
functions permit a simple independent proof of this fact, which we 
shall reproduce here. 

Let 


VP Pu Pe) = =O (m— WY, 


so that 
aV _ 1 ee 
LY =D (,— 9h = D4. 


We shall show that }~ f,<0 in a small region on the sphere }> pi=38 
about (1,1, 1), (except for this point itself), and this will demonstrate 
the lemma. Of course at (1,1,1), }> f,=0. Using the Taylor series 
expansion, we obtain 


(8) DA=D@e-vEw% 
j T Op, 
i = Ra Ff, 4/2 
+ FE he, - ve. - vp ee +00, 


where all of the partial derivatives are evaluated at the point (1, 1, 1). 
Let us simplify some of the terms. The Walras Law reads )> p,f,=0, 
so that differentiating with respect to p,, we obtain 


fi = 
(9) f+ Yeeyt 0. 


At (1,1, 1) this implies }>, @f,/@p, = 0, and therefore the linear terms 
of the Taylor series vanish. Differentiating (9) once more with respect 
to p,, we obtain, at the point (1, 1,1), 
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ye of, of. 
1¢ as ee fee 
(10) T Opp, * Op, ) 


From (7), (8), and (10) we see that 


Dh = —FE @— Were — Yew + 00), 


where 
2C, —C, —C 
(cn) =| —C, 2C, —C 
—C, —C, 2C 
Therefore 


rh= -L tp, ~ p,)' + (p, — Ds)" + (Ps — DE + CV). 


But I claim that on the surface }> pj = 3, with all p, > 0, we always 
have 


Q1) P+ (DI + (PY ESO OH-I. 
Accepting this as correct for the moment, we see that 


DA< —$cv + OV), 


and therefore }>f, is negative for V sufficiently small. 
Equation (11) is demonstrated as follows. First of all, on this sphere 


38>) 0p. 
Multiplying by >> p, and expanding, we obtain 
8 m= Pit pit pit App. + DP. + PP)» 
and by subtracting 9 from each side we obtain 


3(-L pi + 20 p, — 3) > —2(p} + pi + pi) + Ap.p, + Pp. + PDs) 


which is the same as (11). This finishes the proof of Lemma 4. 

We shall shortly exhibit some excess demand functions, derivable 
from a utility function, which satisfy the hypothesis of both Lemmas 
2 and 4. They will provide us with examples of a unique equilibrium 
point, which is unstable locally. Does this imply global instability? 
In a certain sense this becomes a matter of definition, the problem 
being that it is possible that one of the prices might become zero. 
In this event, some care should be taken, since the differential equa- 
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tions (2) no longer describe the motion of the prices. In our examples, 
we shall use special techniques to demonstrate that the paths actually 
stay away from the boundary of the positive orthant. 

Now let us consider utility functions for the first consumer of the 
form 


+1 +1 +H 
U(2,, 2, %) = -(#— . a + eS) , 
2 


xt xy 


where (q@,, @,, @;) is an arbitrary non-negative vector, and a is a positive 
constant. Let the initial endowment of the first consumer. be repre- 
sented by (J,, I,, I,). As before, we permute the utility function cycli- 
cally for the second and third consumers, and also permute the initiai 
holdings. 

We shall show that instability arises whenever a > 1, (a, a@,, a) is 
close to (b, 1, 0) with b>(a+1)/(a—1), and (J, J,, J,) is close to (1, 0, 0). 
(The specific meaning of ‘‘close’’ will be clarified.) This will give us ex- 
amples with none of the disagreeable features of the simple examples 
in Section 2. Let us first examine in detail the case 


(@,, &, @) = (6,1, 0), 
(1, Ly, 1) = (1, 0, 9), 
with 
a+1 4 
a—1l 





b> 


A routine caleulation shows us that the excess demand functions 
for the first individual are given by 


/+a) 
bpt 








= bp/are + pa/a+a) oe 
os P: 
(12) %,= py ara) bpe/it+a) + py ata) ° 
g=@. 


The reader should notice that if '=1, and a— oo, these excess 
demand functions tend to those of Section 2. 
Let us show that the conditions of Lemma 2 are verified. We have 


A= 2:Pv Pv Pd) <Q, 
Op, 


B= O2(Ds, Pi» Pa) >0. 
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Therefore the system of market demand functions (4) has a unique 
equilibrium point given by (1,1,1). In order to show that it is un- 
stable we compute 

"hy —_ of; 

Op, 

at (1,1,1), or 

_ ab—(b+1)-a 

(6+ 1)(a + 1) 





which is positive if b > (a + 1)/(a@ —1). It follows from Lemma 4 that 
for this relationship between a and b the unique equilibrium point is 
also locally unstable. 

The next step will be to show that for this example the prices stay 
away from the boundary of the positive orthant. This will be a con- 
sequence of the following lemma. 


LEMMA 5. Let 


9 = Min (4, Ps Bi), 
Pi Ps Ps/ 


There is a small positive constant K, such that if 


9 <K, then 2 >0. 
dt 
Let us assume that the minimizing value above is p,/p,=c, say. 
The other cases are handled by a cyclic permutation. We want to 
show that d(p,/p,)/dt > 0 or p,f, — pf; > 0, if e is sufficiently small. 
If f, and f, are computed according to (4), we see that p,f,— pf; 
has the sign of 





Ps) 44 Ps Ps 
(13) Pry a) aed Pr. a) < Pr ° 
(2) +1 (2) + ~ Hey +1 
1 a 3 


The first term of (13) is equal to 


(z) +1 
4)" + 1 ’ 





and since p,/p, < l/c, we have 
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Ps 7% 
Ps > c 
a/ (+a) 1 a/(1+a) 
b( Pe te 4 1 
( Ds ) ( c ¥ 
Also we have 
a 
Pi 1 
aj (+a) - Se 
b( Be +1 6 
(2) 


Adding these terms together we see that (13) will be positive if 


(Ay r1-2-Way"" 30, 


and this is certainly correct for small c. This demonstrates Lemma 5. 

These observations, taken together, show us that the system of 
demand functions defined by (12) gives rise to global instability of the 
price adjustment mechanism. 


h, 


3 








R 


FIGURE 2 
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One more remark is in order. I claim that if we take an excess 
demand which is close to (12) on a large proper subset of the positive 
orthant (and which has derivatives which are close) and form the 
market demand functions using the above method of cyclic permuta- 
tion, then this will also give rise to global instability. (Such an ex- 
ample may be obtained by taking (a,, a@,,a@,) close to (b,1,0) and 
(I, L, J,) close to (1, 0,0), or in other ways.) 

In order to see this let us interpret Lemma 5 geometrically. 

Figure 2 represents the surface of the sphere )> pj = 8, and the 
shaded region the set of points with 


Min(2,»,%)>xK, 
P, DP. Ds 

for a small value of K. Lemma 5 tells us that a path which begins 
in this region, with K sufficiently small, will never leave the region. 
This is true because on the boundary of this region expressions such 
as p,f, — p,f,; are strictly positive. But this later fact will also be 
true for market demand functions based on individual demand functions 
close to (12). 

It will also be true that inside this region there will be no equi- 
librium points other than (1,1,1). For if the new demand functions 
and their derivatives are close to the corresponding quantities for (12) 
inside the shaded region, then A for the new demand functions will 
be negative, and since it is homogeneous of degree —1 (being a de- 
rivative), A will be everywhere negative in the region defined by 


m>K > 2.>k. 

Py DP, Ps 
A similar statement with a reversed inequality holds for B, and there- 
fore Lemma 3 may be applie] to show that there are no equilibrium 
points other than (1,1, 1) in the shaded region. 

A similar argument convinces us of the local instability of the equi- 
librium point (1,1,1). It should be remarked that even though we 
have said nothing about possible equilibrium points outside the shaded 
region for the new demand functions, a path beginning inside the 
shaded region will not leave this region. 
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SIMULTANEOUS LEAST SQUARES: A DISTRIBUTION 
FREE METHOD OF EQUATION SYSTEM 
STRUCTURE ESTIMATION* 


By T. M. BRowN 


SUMMARY 


A METHOD of structure estimation is here developed which separates 
the regular from the random component of structure by minimizing 
the vector distance between the observed data and the corresponding 
values computed from the systematic structure of a complete model 
treated as a whole. The method is, accordingly, the simultaneous 
equation counterpart of single equation least squares. It too is dis- 
tribution free. 

The sample error variance-covariance matrix is constructed. Its 
elements are zero for an infinitely large sample. This, plus the 
hypothesis that only true structure could produce the above minimum 
vector distance, suggests consistency. The method is tested on a 
model with known structure, which generated its own data and 
moments. The estimates are equal to the true values. 

The method is also tested on a model with unknown structure, 
using actual data. The results are compared with estimates by other 
well-known methods. Further research on the relative efficiency of 
estimation methods is suggested. 


1. INTRODUCTION 


There are three uncomfortable assumptions involved in the method 
of full information maximum likelihood (FML) for the estimation of 
economic structure. One is that the disturbance vector follows a 
multivariate normal distribution. A second is that a large sample 
of basic data is available for estimation purposes. And the third is 
that the best basic principle on which to base the estimation procedure 
is that of maximum likelihood. 

With regard to the first assumption, while the distributions are 
probably bell-shaped and may tend asymptotically to normality, since 
each disturbance tends to combine a large number of lesser random 
influences, the normality assumption may not be closely met at all in 


* Manuscript received February 26, 1960. 
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small samples. But small samples are the rule with economic time 
series appropriate for structure estimation of complete models. The 
passage of time may not correct this, for structure evolution may 
force the gradual rejection of data from periods too far in the past. 
This tends to invalidate the first two assumptions. 

With regard to the third assumption, while optimal properties have 
been proved for the maximum likelihood (ML) approach to structure 
estimation, optimal properties have also been proved for the least 
squares (LS) approach. This produces a conflict of loyalties, a dichot- 
omy. To demonstrate this point, let us assume that the disturbance 
of a single linear relation satisfying conditions of the Markoff theorem 
on least squares (David and Neyman [6]) is normally distributed. 
Then the ML and LS approaches lead to identical results, which are 
best linear unbiased estimates. There is no conflict. But let us next 
assume that the disturbance follows a non-normal, say a Gamma, 
distribution. Then the ML and LS estimation equations are quite 
different. Yet the LS estimates retain the property of -eing the 
best linear unbiased estimates, and in addition are much easier to 
compute. In this case which method, ML or LS, will produce better 
estimates, if we ignore computation costs? 

In the case where the parent distribution of the random variable is 
not known, there is no conflict. LS must be used. It can be used 
with equanimity because it produces best unbiased estimates of struc- 
ture. In addition, its sampling error formulas are appropriate for 
samples of any size. 

The LS method appeals to the principle that minimizing the vector 
distance between the observed data and the systematic component in 
its explanation produces optimal structure estimates. Intuitively, one 
might assume that this is a necessary condition on the true structural 
values, when the universe of data is available. Some regular com- 
ponent added to the random component because of imperfect estima- 
tion of structure would increase this vector distance’. Conversely, it 
would be impossible to decrease this vector distance between true 
systematic and observed data by shifting some of the random com- 
ponent into the systematic. 

These same concepts can be readily carried over to the complete 
system of simultaneous equations comprising an economic model. Let 
Y, Z, and u be vectors of endogenous, predetermined, and disturbance 


' Let we be the true random component and let r; be part of the regular structure 
which is added to yu; because of faulty structure estimation. Ely:+r:)?=Ey:+2Euere+ Eri. 
Let d= 3 (v2? + 72), ®@ = Dw; d*> dG. Note that Eu:r:=0 since r; is part of regular 
structure. 
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variables respectively, in a complete economic model. Let (YZ) = X 
and « = X— X. Let B,C, and S be the matrices of unknown struc- 
tural parameters in the model which are to be estimated, such that 
(1) By' + Cz =w; 


T 
E[u'u] =>; Lowi = §; Av’=wu’; A=[BC]; B is square. 


This complete system generates the endogenous variables from the 
structure, the predetermined data, and the disturbances through the 
process 


(2) y’ = —B'C2' + Bow’ = Fe’ + u,'. 


This is the simultaneous solution of (1) and produces the final values 
of y after all the interactions among the y have been taken into 
account. 

The complete explanation of y separates into a regular component 
y-., computed from the regular part of the structure, plus a random 
component wu, arising from the random disturbances in the system. 


Thus 


(3) y,= —B "C2 = Fr, 
(4) “un Boy’ , 
(5) Y¥— Ye = U,. 


Equations (2) and (8) may be referred to as the structural reduced 
form of the system, the reduced form based on the structural parame- 
ters. Since y and z represent observed data, uw, represents the dif- 
ference between the observed data and the component explained by 
the regular structure for the economic model treated as a whole. 
Parallel with the least squares principle for a single equation the fol- 
lowing questions are pertinent: (a) What values of B and C will cause 
the vector distance of u,,(t = 1,2, +++, T) to be a minimum? (b) Will 
such values of B and C have optimal properties, as in the single 
equation case, e.g., will they be consistent? Question (a) can be an- 
swered by a well defined mathematical process which produces unique 
results’, With regard to (b), using the same kind of intuitive reason- 


2 The minimization of a vector distance in the space of endogenous variables suggests, 
in analogy with the theory of orthogonal regression, a lack of uniqueness under transforma- 
tions of the endogenous variables. For any given units of measurement the results are 
unique. It can be shown that the results are invariant under orthogonal transformation 
of endogenous variables but not necessarily under any transformations in finite samples. 
In just identified models the results are invariant under linear transformations generally 
I am indebted to Professor H. Theil for pois.ting out the problems of invariance under 
transformation. 
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ing as that above for the single equation case, we would expect to 
find consistency.’ 

Because of its parallelism in concept to single equation LS, the 
method suggested can be referred to as Simultaneous Least Squares 
(SLS). It takes fully into account the simultaneous interactions of 
all of the endogenous variables in the complete system. It obeys 
‘‘Markoff conditions’’ within the multi-equation system, since it works 
with the y’s as linear functions of z’s and w,’s only, where the z’s 
and u,’s are statistically independent. 

The vector distance which is to be minimized is d, where 


T 9 
(6) @ =>) >) tis « 

t=1 t=1 
There are g encogenous variables and equations in the system, k pre- 
determined variables, and 7 time periods, for each of which ¥y,, z, 
are observed. We also have 


T 
(7) d= D=tr)) uiu.,, = trTS, = trS} . 
Thus the SLS method involves minimizing the trace of the variance- 
covariance matrix of the structural reduced form disturbances. 

This reveals that SLS is closely related to FML, or to Quasi Full 
Maximum Likelihood (QFML), in the case where the FML procedure 
is used when the distribution of the w’s is unknown. For QFML in- 
volves finding the values of B and C (in the reduced model with all 
identities eliminated) which maximize the function 
(8) La) = lox 5 
({11] pp. 212-3; also [3]). These values will also maximize |B|?/|S}, 
and hence will minimize |B|-|S|!B|-. By the multiplication theorem 
for determinants |B | B-'| = 1, and |B-||S||B-”| = |B“"SB-” | =det S,. 
Thus QFML minimizes det S, while SLS minimizes tr S,. 

Since det S, includes all of the covariances as well as the variances 
in S,, while trS, includes only variances, it might be thought that 
this would in itself make QFML a more comprehensive method, using 
more information in the estimation procedure. But note that trS, = 
tr B-'SB-* involves all of the covariances as well as the variances of 
S, and hence includes all of the covariances of the primary disturbance 


* Consistency is the property that the estimates converge with probability 1 to the 
true values, as the sample size becomes infinite. A formal proof is given in the accom- 
panying note in this issue by M. Nakamura. 
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vector wu. Also min det S, permits large variances to be offset by large 
negative covariances occurring in the right positions. Min tr S, contains 
only the (positive) variances in S, and, in addition, appears to have 
more meaning than min det S, in relation to the estimation problem, 
as discussed above. 

Which of the two methods is in fact superior can only be decided 
by further research. Questions of consistency (QFML was proved te 
be consistent in [11] Chap. II), sampling errors, and computational 
cost must be studied and weighed. 


2. THE SLS ESTIMATION PROCEDURE 


In the QFML method (cf [3], [4], [11]) it is necessary to eliminate 

all identities contained in (1). For SLS this is not necessary, but it 
seems advisable to do so anyhow, to reduce the size of the matrices 
in the computations. The reduced model with identities eliminated 
can be represented by 
(9) Aja =u, or By +C,2z.=u,. 
Assume that there are g, elements in y, and k, elements in z,. Let 
the unrestricted or free structural parameters in A be (a) = (a, +++, a,) 
and let A also contain the known or restricted parameters 1, G,, @,, 
and a, = ka,, where a, is contained in vector (a). 4,,da,, and k are 
assumed known in advance, possibly from cross-section and engineer- 
ing measurements. In the estimation procedure they are always used 
with their known values. The individual elements a,, of A, are func- 
tions of the parameters of A. 


(10) a, = g(a, 9°, Q,; 1, As, a, k) = f(a); (é = he cere m) . 
From (9), we have 
vv, + Bo'C,2}.= Bru} = uj, . 
Then 
? 
(12) D,=trS, =tr Dy — F.zi\y, — z,F,) = Da,) = Da). 


We now drop the subscript: r for the sake of economy, but it is 
understood throughout. We also adopt the notation for moments. 


(18) m,.,=Soaes % = X,— Xs [m,.,) = Mos 

{i,j =1,2,+++,(9+ b}. 
Then (12) becomes 
(14) D(a) = tr(M,, — M,,F’ — FM,, + FM, F"’) . 
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This is a curved surface in the space of (a,, -++,a@,; D) where a,, +++, a, 
are the unrestricted free parameters of A in (1). 

Our problem now is to find the values (@) which bring D(a) to its 
minimum value D. An iterative process is used. An initial vector 
(a), somewhere near (@) is obtained from some simpler single equation 
method. The direction of steepest descent from (a), is given by the 
vector —(@D/da),. 

To obtain fast convergence (or even any convergence) this vector 
direction must be corrected from steepest to shortest descent. This is 
done by bringing it closer to a geodesic line through (a), and leading 
toward (a). This correction is achieved with the matrix [@°D/6a’da);*=D, 
(Cf. Brown [3], Chernoff and Divinsky [4], and Crockett and Chernoff 
[5]). The iteration formula is 

1 
(15) (a), = (a) + Aa); (a), = —W(2?)/_22 |" 

As in the FML case, it is to be hoped that some simple component 
matrix, say P of D,, will be found to act as a substitute for D,, 
providing good convergence at reduced cost. 

2.1. Formula for @D/éa. Formula (14) suggests that we first find 
dD = (aD/af)df’, where f = vec F. From (14) 


dD = tr {(—(dF)M,, — M,,dF’ + FM,,F’ + dFM,,F"} 





(16) = 2 vec (FM,, — M,,) vec’ dF , 
(17) 4 = 2vec (FM,, — M,,) . 
Then 
(18) aD _ 4 

da of _ 6a 


As D is approached, @D/da—0. Note that from (17), (3), and (5), 
using z as a vector of instrumental variables, we obtain, 


(19) 9D — 2 vec Soule. 


af 


This approaches zero for a large sample, but need not be zero for 
a small sample. For a small sample SLS does not force @D/df to be 
zero, and @D/@a is zero by virtue of the product of the matrices in 
(18). But note that in the case of the statistical reduced form the LS 
regression of each y on all z’s, FM,, — M,,, is equal to zero for samples 
of all sizes. This source of small sample error in all single equation 
LS estimates may not be present in SLS. 
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2.2. Computation Formulas for [df /da’}. 
(20) F, —= — BC, ’ 
(21 dF, = —B;'dC, — B7'dB,B>"C, , 
(22) oF, _ _p- { r ac, 90, B;| 9B. jp, 
da, da, L 0a, 





where a, is an element of (a), not an element of (a,). 

aC,/da, and @B,/da, are easily computed from the structure of the 
reduced model. Numerical values of B;' and F, are readily available 
with each iteration. With these f/@a can be computed for each itera- 
tion. When working with a specific model, however, much time will 
be saved by applying (22) to the model structure and working out the 
algebraic formulas for each 2f,,/0a,. 

2.3. Computation Formulas for D,. The direct differential approach 
to this matrix is tedious and impractical for large systems. However, 
a partial differentiation approach produces matrix formulas which can 
handle systems of any size. 

Suppose x = x(u, v), y = y(u, v), then 


y ( 


y 
4 ou 


(23) (ev) = «(2 Ju . 


The case is similar for 0(xy)/dv. 





(24) dD =(2 7 Pas" df’ = [ aa’ dD = PLE |e. 
(25) we elles: 


Applying (23) to (25) we have 




















to, or) a 123 ae Ye wi Ea 
(26) 
L(PVED- BSG PIE 
20 (Sp) = apap ht + apap 
{ = [ |X lao , 
df Of JL da J 


Hence, 








180 T. M. BROWN 


° £Q2) - [S814]. 





From (25), (26), and (28), 
@) FA a (Fee | . lee lZ] aS +2. 


In (29), |[@f/@a’da]| is a three dimensional array with gk rows (of 
the order of f), » columns (of the order of a), and n layers (of the 
order a). It is essentially a three-way matrix and could be designated 


as a tri-mat. The product of the vector @D/@f and the tri-mat is a 
matrix Q of order » by n. It is computed in row by row sequence 











as follows: 
Row 1 of Q = 2D or] 
Of _da,da 
(0) : : 
Row n of Q = 2P2f or. 
df _da,0a 





The unknown portions of (29) so far are @f'/da’da and #D/df'df. 
From (21) we know that @f/0a consists of elements like b”f,, from B?* 
and F. From these terms we build up the elements of 0°f/éa’da, 
using formula (23). Terms like 0b”/da,, are computed from 





(31) dB = —B-dBpp-, 2B” _ —B-| 2B lp, 
00, 04., 
@D/af'af being easily derived from (17). 
(32) or = 2vec FI, & M,, — 2 vec M,,, 
#D 
—— = 2], M,, . 
(33) af oF 7® M,, 


2.4. Pas a Substitute for D, for Iterations. The P component 
of D, has all of its parts readily available and hence suggests itself 
as a convenient substitute for D, in the iterations. In the tests on 
Models X and Y below it proved to be satisfactory and gave rapid 
convergence. 


8. THE SAMPLE ERROR VARIANCE-COVARIANCE 
MATRIX S(@) OF THE SLS ESTIMATES 


This section is sketched very briefly, to conserve space. Let E be 
the expectations operator, and d, be the sampling differential operator, 
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as we assume theoretical repeated samples of Y and an estimate (a) 
corresponding to each sample. Z data are assumed constant in the 
repeated samples. 


(34) da=a-—  E(a). 


Consider the change in D as w’ ange from sample to sample. Omit- 
ting the bar and subscript s 


(5) 4D = Fda, + oT oni 











ll Yor 
aD\ #D aD #D 
AS) oe Sh ae, + oo ee, + 
(36) Ge) ba,0a, * Gaga," * tae” 
free OO dy,,=0, as D>D,a—a@; 
a,0Y yr 
2 | a, ’ | @’D |a ‘=0 = eee eee eee 
(37) FS aba a + 6a'dy; . s¥1 (Yr Yury Yirs You »Yor)s 
aD aD 
'da} = D;| | pe 
(38) E\de'de) @a’d dy;da 
) a ‘ee 2 a 
(39) E{dy;dyr) =| . ° , 
Pete fie _— 


which reduces to I,@ >>, when we assume zero serial- and auto-cor- 
relation in the disturbances. From (17) and (18), 


aD _w Vo’ 
(40) <0 = 2vee (FM, M,)| af ] 


, 


Yu *** Yn @u °°? 2x | 


(41) a,(57) ojea. 8. 





ere: nae. 


Yir °°* Yor Zr *** Zer 








&’D OY+ 7 Of’ of’ , 
= — 2 vec —*Z,— = — 2NI,® Z,~— , 
(42) ayia oy, "8 0&8 om 
‘@D ie oY, , of’ 
dyu0a Svee dyn" 0 da 
#D =-—2 vec a¥ry af” 





dyorda yor = aS 








182 T. M. BROWN 


where N is gT x gT.The small sample formula is 
(43) S(a) = 4D; 21, @ Z;N'ly @S,NI,@ ZLD> . 
da’ da 
For large samples we change the moment formula (13) to 
T 
>> U4, « 


S(a) now becomes dimensionally of the order (4/T)V, showing that 
S(a)—0 as T— 0, 


(44) M,,2, = 


4. GOODNESS OF FIT 


All of the conventional goodness of fit statistics can be calculated, 
once the model has been estimated. A suggested selection is: t, = 
a,/S(a,); coefficients of variation, V,= {S(u,)/Y,} x 100 and, V,, = 
{S(u,,)/ Y,} x 100; coefficients of correlation, R,=[{S*( Y,)—S*(u,}/S(Y,)f”, 
and R,, = [{S*(Y,) — S*(u,,)}/S*(Y,)}’; and, finally, von Neumann ratios 
for each u, and u,,. 

It is also of interest to consider the possibility of a single measure, 
corresponding to R, of the goodness of fit of the complete model treated 
as a whole. Recent work on this problem has been done by Hooper 
[8] and Theil [13]. In terms of the SLS estimates, the measure 


(45) R, = (= Ti =.) 


gives a parallel to the R,, above. This formula, however, omits covari- 
ances. This omission can be remedied, if desired, by using determi- 
nants rather than traces in (45). We recall that trS, already includes 
all of the covariances of S. An interesting variation which also in- 
cludes covariances follows from Hooper’s formulas and his suggestion 
on p. 254 of [8]. 





(46) R, = {Su a“ M3s,|* 


5. COMPUTATIONAL EXPERIMENTS 


The SLS method was first tested on Model X, the true structure 
of which was known. This same model also had been used for test- 
ing in [3]. In this model @ and J are given. Also assumed given is 
a fixed sample of Z;. From these Z, data, Model X generates its 
possible universe of Y data, from which population moments could be 
computed. From these moments experimental computations were made 





SIMULTANEOUS LEAST SQUARES 183 


to estimate the structure a and 3, using the SLS method. If the 
derived estimates equal a and , this would demonstrate the consistency 
of the SLS method. At the same time, the number of iterations to 
convergence and tic general volume of computation could be appraised. 


5.1. Model X. 
(1) C= .5+ 8W+ 4r+4u,, 


a & a; 
(47) (2) W = , + oY + .) + Us ; 
0 1 2 


(3) W+xr+T7,=Y, 
(4) C+E=¥Y. 


Endogenous Variables Predetermined Variables 
C = Y, = consumer expenditures T,= Z,=government net revenue 
W = Y, = wage-salary bill E = Z, = all non-consumer spend- 
xz = Y, = non-wage income ing on newly produced 
Y= Y, = total income, production, final goods 
and expenditure Y_, = Z, = value of Y lagged one 
time period. 


With identities eliminated, and variables converted to deviation from 
mean values, the reduced model becomes, 


(®) &5= i iit fh oie ”} C, = i , sf BC, = A,; 
(6, — 1) b, b, 0 b 
& = (Yar Vas Zi» 2 %) = (X — X); w= (—w, —u,); Aw’ =u’. 


We assume the universe of u to produce a constant variance-covari- 
ance matrix. 


(49) > = Elu'u) = free pal 
.004 .012 


To make this matrix more comparable to small sample formulas, we 


multiply it by a scalar 7 = 25. 


(50) >* = TE=|" ‘al: 

et 3 
Moment formulas are defined as follows: 
(51) M:.2, = E(X, — X)(X, = X;) 


is a population moment. 








184 T. M. BROWN 


m,.,= 2 (X,— XX, — X) 


is a small sample moment, with 7 observations in the sample. But 
for this experiment we also use M2, = Tm?» with 7 a scalar equal 
to 25. The z moments as derived from the Z, data are assumed to be 


7 a ee 
(52) M,, = (m,.,) = TM?, =| .1 2.0 1 
2 .1 8.0 


The moments for the universe of detailed Y data which the above 
structure would generate from the given Z, data can be computed 
directly from the formulas: 


(58a) M,, = TM}, = B>">>* B>” + B7'C,M,,CiB;” , 
(53b) M,, = —B;C,M,, . 


The subscript r is henceforth usually dropped from A,B, and C, 
but is understood. Also the lower triangle of all symmetric matrices 
is omitted. From formulas (53) we obtain 


(54a) M,, = a as 2.159375] 
Ys 6.708819 
zy 23 23 
(54b) M.. = ¥ —-187500 1.250000 pool 
" — y| —1.454167 2.750000 —.420833 


The free or unrestricted parameters of Model X are evident from 
(47). They are 


(55) (a) = (a,, a,, 5, b,) , 
(a) = (.8, .4, .8, .2). 


The same opening vector point from which iterations began in the 
FML experiment [3] is used in this experiment with SLS. Thus 


(56) (a), = (.6, .3, .4, .1), 


with percentage deviations from true values (—25.0, —25.0, +33.3, 
—50.0). These deviations are of the order of 334 per cent. 

The caleulations for Model X were done by Mr. J. A. Auréle Leduc, 
who took considerable interest in the work and made helpful sugges- 
tions. In the main, the calculations were carried out to ten significant 
figures to keep rounding errors from accumulating into the first four 
or five significant figures. The data were scaled so that ten significant 
figures usually meant nine decimal points. The numerical results 
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which follow are, however, shown only to six decimal points to save 
space. Matrix inversions were done by either adjugate or triangular 
methods, with the square root procedure applied in the case of sym- 
metric matrices. Sum checks were used throughout for all matrix 
products and inversions. 

i. The First Iteration: 


(57) 7 ee Poe 700000}, a =a 1 07. 
mess ° = : 
.600000 —.400000 0.1 
(58a) pen an - 689655 geod, 
1.034483 —.689655 


(58b) —B=C,, = F, = | —.206897 .689655 - reat 


—1.310345 1.034483 —.068966 


69) H=FM -( —.113798 1.370690 pond 
ee na 


—1.220690 1.931034 —.365517 


(60) tie (2) = — vec L, = vec (M,, — H,) 
2 \ Of /o 
== (— .023707, -—.120690, .372845, -—.233477, .818966, — .055316). 


(a) Formula for |df/da’|. Applying formula (22) to Model X gives 


bE 
|] os b* {(1 + fut fu)s (fu + fax), (Sis + Fua)} | 
da 
61) 0 0 b? 
b"F 
le {( + fu + fu), (Fis + Fa); (fis J Fx)} F 


0 0 b” 


where the 5” are the elements of —B;' and the f,, are elements of F’. 


ow [2] 


= [ —.142687 _—-«.475624 .083234 —.21403]1 713436 .124851 
— .903686 .713436 —.047562 -—1.355529 1.070155 — .071344 

— .624257 2.080856 -062426 .356718 —1.189061 — .035672 

0 0 1.206897 0 0 — .689655 
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(b) Vector of Steepest Descent. 
qa -1. (82) = -2 (DY) 


2 \@a 2 \af/loa 
= (.604358, 1.114437, -—1.268175, .488134), 
: | 2 = 1,QM, - 
(63) 2 Lof’af 1@ M,, 


(c) The P Matriz. 


+> "(Sl dee hleh 

2 da’ 2Lafafs Loa 
(64) 1.570816 2.380793 .290401 . 046382} 
5.599859 .468710 —.043462 
11.689263 .483085 
5.796671 
Segoe —.762854 —.013334 — .010082 
.508161 —.001594 .009989 
.086179 —.006344 
.178214 





(65) 2P;° 








(d) Vector of Shortest Descent. 
(66) (2?) Pe = d, = (.243125, .106601, —.122221, .092197), 
a /o 
(67) (a), = (a), + hd, , 
with h, selected at 0.5 = (.721568, .353301, .338890, .146098). 
(e) Distance Formula for Checking Progress of Iterations. From 
(14) we have 
(68) D = trM,, + vec (FM,, — M,, — M,,) vec’ F 
= trM,, + vec (L — M,,) vec’ F, 
(69) D, = 1.527716; D, = D = 1.059028 . 


Progress in the iterations can also be checked from the changes in 
aD/af and aD/da. 

ii. The Subsequent Iterations. To save labor, the P, matrix was 
used for three iterations, and then a new matrix P, was computed 
for subsequent iterations. As the iterations proceeded, Mr. Leduc 
suggested that @f/da’ also be kept constant for three iterations, to 
shorten the computations. In testing this proposal it was found that 
the temporarily constant df/da’ had the further advantage of giving 
faster convergence. 
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After six iterations, using only two P and two @f/da’ matrices, the 
estimates (a), had converged to the true values to at least six decimal 
points. The sequence of h values selected for the six iterations were 
(.5, .5, .75, 1.0, 1.0, 1.0). 

Thus, we have, at the minimum, an empirical test that the SLS 
method, in minimizing the vector distance (y, — y,,), correctly separates 
and estimates the systematic and random structures, given the universe 
of data. Also, it was noted that at D, 4(@D/0f) = tr(FM,, — M,,) = 0, 
a necessary condition on the true structure and population moments, 
(cf. (19)). 

iii. Truncated Simultaneous Least Squares (TSLS)’. As in the 
case of FML [3}, the cost of the SLS procedure can be considerably 
reduced by using a truncation procedure. This gain will, of course, 
be at the expense of a certain bias or loss in accuracy. The trunca- 
tion tested for Models X and Y was to use a constant [0f/da], and P, 
matrix for four (or three) iterations and to terminate computations 
at this point. The results for four iterations on model X are shown 
in Table 1. 














TABLE 1 
TEST OF MODEL X STRUCTURE ESTIMATION BY TSLS 

True Values Estimates by TSLS e-TSLS 
a 8 .802339 +0.3% 
a 4 398539 -0.4 
a3 . 298840 -0.4 
a m . 196402 -1.8 
D_ 1.059028 1.059137 +0.01% 





The bias or error formula used in Table 1 is 


(70) e, = 4%—% x 100%. 
a, 

iv. Test of Professor Theil’s Two-Stage Least-Squares (TL) Method. 
Before leaving Model X, it seemed like a good opportunity to test 
Professor Theil’s interesting new single equation method [13]. It gave 
perfect estimates on Model X, verifying its consistency, for which 
there is, of course, theoretical proof. (See also Basmann []1)). 


5.2. Model Y. 
Model Y is a small model of the Canadian Economy, built from 


6 The initials TSLS are not here to be confused with those for “‘two-stage-least-squares,”’ 
dealt with in the next section. 
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actual economic data which are now considerably out of date. The 
data used are for the years 1927-41 and 1946-52, comprising 22 ob- 
servations. Though out of date, Model Y is useful for continued ex- 
perimentation, since it is small and based on a sample of actual data. 
Also, since its data and structure remain constant, it permits com- 
parison of estimation methods. It was used in this way in the ex- 
perimentation on FML reported on in [3]. It will be recognized that 
Model Y corresponds, except for one variable, with the model reported 
on in [2]. This latter model covered only the years 1927-41 and 
1946-49. 

In this article Model Y is described in detail and its moments are 
given. This is done to enable readers to check the results shown 
here and in [3], and to-perform any further testing which may occur 
to them. 


MODEL Y 
(1) C= b,, | a + bs Y. + ¢,C_, + Ciel -1 + Cy th, 
(71) (2) Y,=6,(Y + S,) + ¢,(Y + S,)., + Cus + Oo + tr, 


(3) Y=Y,+Y,, 
(4) Y+S,=C+S8S,+S,. 
All variables are measured in billions of 1935-39 constant dollars, 


except for t, which is measured in years (calendar year minus 1926)’. 
The variables are defined very briefly as follows: 


Y, = C = consumer spending on goods and services, 

Y, = Y,, = disposable wage-salary income, 

Y, = Y, = disposable non-wage or property-enterprise income, 
Y,= Y+S,=GNP, 


Z, = C_, = C lagged one year, 
(72) Z, = L,.-, = liquid asset holdings of households at end of pre- 
vious year, 
Z; = (Y + S)-. 
Z, = t/10, 
Z,=S,= Y-—C, 
Z, = S, = GNP — Y. 


The moment matrices, defined as in (13), are 


® Mr. Lou Scott added his computational abilities to the work on Model Y, and his 
neatness, accuracy, and helpful co-operation are appreciated in enabling the work to be 
completed in the time available. 
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el Y2 Ys Vs 
y:{48.017285 38.077358 21.738878 77.633277 
M,, =" 30.494605 17.174202 62.191638 
Ys 10.933463 35.894292 |’ 
Ti 128 .449409. 
2 22 Z3 zs % % 
yif 45.083213 66.468980 74.176188 23.989559 1.798951 17.8170417 
M... = ¥2| 35-676996 52.058769 59.058006 19.386505 9.591449 14.522831 
™* ys} 20.163852 29.848450 33.532672 10.367895 6.368787 7.786627 |’ 
(73) ys72.553932 103.353066 120.752312 38.814432 20.452653 30.363479_] 
21 a Z3 2, % r2,) 
21 [ 43.196648  62.084527 69.908005 22.580436 10.757635 ple eed 
Za 104.050177 101.694914 33.713850 15.438241 21.445845 
u.=-* 116.458406 36.314400 18.414490 28.161634 
a= 13.818636 5.764841 9.060032 |" 
Zs 4.161285 4.492417 
tol. 8.054021 | 








The results of some of the various methods which have been used 
for estimating the structure of Model Y are listed in Table 2. Some 
methods not here listed are reported in [3]. SLS results are put in 
the lead position and error formula (70) is adapted by replacing a, 
by a,,, the SLS estimates. Either these or the QFML results are as- 
sumed to be the best available estimates of Model Y structure. In 
the computations on Model Y constant A values of 1 were used through- 
out, and six iterations gave 0D/da=0 up to six decimal points. 
Iterations were started from (a), = LS estimates, just as for the QFML 
estimates in [3]. P, and [@f/da’], were used without change for itera- 
tions 1, 2, and 3, and then P, and [@f/0a’], were used for iterations 
4, 5, and 6. In computing TSLS as a short-cut method, three itera- 
tions with the original P, and [@f/@a’], were found to give a closer 
approximation to the SLS values than four iterations. The results 
recorded in Table 2 suggest that if optimal single equation results 
were used for (a),—for example TL or LI (Limited Information)—TSLS 
would give results very close to those of SLS. 

In this small sample case, SLS converged as rapidly and easily as 
for Model X. The vector @D/df, however, did not reduce to zero, 
although most of its elements became small. On reflection we realize 
that (19) should not usually be zero for a small sample. Our method 
does not force this equality, and this may be one of its strengths. 
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TABLE 2 
COMPARATIVE ESTIMATES OF MODEL Y STRUCTURE 
Structural! 15.6 /QFML-10| ¢- e- 

phen Iterations | Iterations QFML TSLS |TsLs TL e-TL Ll e-Ll 

ers (a) 
biz .678681 .680868 |+0.3 .684290 |+0.8 .682987 | + 0.6; .686677 | +1.2 
bis .064848 .063090 |—2.7 .063965 |-1.4 .087960 | +35.6| .065352 | +0.8 
Cu .351194 .350868 |—0.1 .346262 |-1.4 .335461 | — 4.5] .341344 | —2.8 
C12 -070813 .070369 |—0.6 .071283 |+-0.7 .071707 | + 1.3] .072836 | +2.9 
bes . 159513 .159723 |+0.1 -161895 |+1.5 .163828 | + 2.7 .161933 | +1.5 
C23 . 240246 .241150 |+0.4 .236802 |~-1.4 .239683 | — 0.2 .241373 | +0.5 
Cu -032488 .032147 |—1.0 -032782 |+-0.9 .031289 | — 3.7| .031377 | ~3.4 






































6. CONCLUSIONS 


The SLS method of structure estimation for econometric models has 
behaved well in the above tests. Its results compare favourably with 
QFML. It also seems to be easier than QFML to compute, both in 
individual iterations and in the number of iterations to convergence. 
Like QFML, it is a full information method. 

SLS has appeal on theoretical grounds, for intuitive theory suggests 
that only the true systematic and random structures can produce the 
minimum vector distance between observed data and the systematic 
component of the explanation of these observed data. In following 
this approach, one need make no assumptions about the distribution of 
the random variables, or about the sample size. The method is accord- 
ingly both ‘‘distribution free’’ and ‘‘sample size free.’’ 

The present study suggests that it would now be useful for research 
to attempt to decide on the superiority of SLS or QFML. This would 
hinge partly on the relative sizes of S(a) for the same model estimated 
by the two methods, and partly on any predisposition of either to 
produce small sample biases. 

It may only be possible to compare these magnitudes when they 
are in numerical form, related to a specific numerical model. To ob- 
tain such numerical tests, synthetic models like Model X could be set 
up, and samples of data could be drawn from them by Monte Carlo 
methods. The synthetic model used could have the same specifications 
as some model in which the researcher had a particular interest. 
Estimates of both a and S(a) could then be made by each of the two 
methods. It should be possible to work out the sample error formula 
for S(a) of QFML, proceeding exactly as above for SLS, and merely 
replacing D=trS, with D = det S,,. 

Any tendency of either method to a small sample bias and the size 
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of such bias as an empirical function of sample size could be estimated 
by drawing repeated Monte Carlo samples in groups of varying sample 
sizes from the model. For each sample the structure a could be esti- 
mated, preparing the way for estimates of the probability distribution 
of a for each sample size. Electronic computers make research of 
this kind a fascinating possibility. That method which reveals the 
smallest small sample biases and the smallest sampling error variances 
and covariances would win the first phase of the contest. 

The second phase would consist of a comparison of computational 
cost. An efficiency function varying inversely with sample bias, sam- 
pling errors, and computational cost could then be set up to help choose 
the appropriate estimation method for a particular purpose. In mak- 
ing such efficiency comparisons it might be appropriate to include also 
the truncated and the single equation methods. 


Royal Military College of Canada 
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A NOTE ON THE CONSISTENCY OF SIMULTANEOUS 
LEAST SQUARES ESTIMATION* 


By MitrsuGU NAKAMURA’ 


VARIOUS PROBLEMS connected with structural estimation through 
simultaneous least squares procedure have been studied by T. M. Brown 
in a preceding article in this issue’. In this paper he suggests that 
simultaneous least squares estimation has the property of consistency, 
as well as a lighter computational burden as compared with ordinary 
full-information quasi-maximum-likelihood estimation, with which his 
method may naturally be compared under fairly general conditions, 
without assuming any particular type of distribution. Brown shows 
the results of several experiments which produce favorable compari- 
sons in plausible situations. Rigorous proof of consistency by this 
procedure, especially concerning sufficient conditions, is, however, open 
to question. In this note I shall try to present a proof for consistency 
of simultaneous least squares estimation starting from fairly general 


assumptions. 
Our structural model is as follows: 
(1) By, + Tz4,=%u, (¢ =1, 2, -+), 
with assumptions: 
(2) |B) #0, 
(3) E(u,|z,) =0, 


where y, and u, are column vectors with m components and z, is a 
column vector with m components. Our model is complete, having m 
endogenous variables (¥;,,Yx,°***,Ym:) and m exogenous variables 
(2: Za, ***s Znt)- 

Introducing new parameters JJ and new random disturbances v,, by 
definition, we have: 


(4) Te-BT, 
(5) v, = B-4u,. 


* Manuscript received, January 15, 1960. 

1 I am much indebted to Professor L. R. Klein and Professor T. W. Anderson for 
reading the manuscript and making useful suggestions. The writer alone, of course, is 
responsible for any errors which may remain. 

3 See pp. 173-192. 
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We now can obtain transformed relations, i.e., the so-called reduced 
forms: 
(6) Y, = lz,+,. 

We assume the following conditions in the above structure: 

i) Structural coefficient parameters, B,/, are identifiable from a 
priori information about B, /’, and the square matrix B is non-singular. 

ii) Conditional distributions of (u,, u,, +++, Uy) given (2, 2, ***, Zy) 
are independent of the given value of (z,, z,, +++, Zy). 

iii) u,(t = 1, 2,-++) are random variables mutually independent with 
respect to different t and have bounded variances of their components. 
iv) |z,|(¢=1,2, +++) are bounded, say |z| < K, and are not zero. 

v) lim|z,| #0. 

vi) Define ‘‘Limit Directions’ of z,(t = 1,2,+++) as a set of limit 
vectors of z,/|z,|(t = 1,2,-++). Then we assume here that z, has at 
least n mutually linearly independent limit vectors in ‘‘Limit Directions’’ 
of 2,. 

We denote the sum of squares divided by the sample numbers as 


N 
>> (y. — Mz,)'(y, — M2.) = SHU» Yar °° *%s Yws 2» °°» Sy M1) 


or more briefly 
= Si(Yy, 1), 
and the true coefficient structures as B,,/,, and //, respectively, or 
],= —B;'T,. We define simultaneous least squares estimates of /] 
from a sample of N, fi yw, a8 follows: 
N N e ™ 

Min {9° (uy. — Mz)(ys — Med} = So (ue — Aye — Heed), 

neo, \ N t= N t= 
where 2, is a set of matrices —B“J in which B and J satisfy the 
a priori restrictions. We define |/7| for any // as 


imj=(Ssm) » 


i=l j=l 
where 
(1,;) = M1 . 

THEOREM: Simultaneous least squares estimates Il, tend to true 
values I], in a probability sense under assumptions (i) through (vi) 
above; that is, for any given positive numbers ¢,r we can find a 
certain integer N, such that 


Pr {| 1y — 1,|<¢} [1-2 
holds for any N= N,. 
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Proor: Define A/J for any J] as 
All = IT = I], . 
We have the following identity: 


NS}(Yy, 1) = Sy. — Me, — Allz,)'(y. — We, — All2,) 


= > (v, — Allz,)'(v, — Allz,) 


t=1 


N . 
= > {viv, — 2viAllz, + (Allz,)'(Allz,)} , 


or 


N N 
7) Si(Yu, M) = S3(Yun Me) — 2 So vialle, + yo(alle, (Alla). 
The second term in the last expression above satisfies relations 


= ualle,| = = 2| 1 DE wAMz, 


t=1 » t=1 


Ss 25° lan Lg vuts| 
N t= 


{=1 


<x |an|| 2g 


{=1 t=1 





< 2|A1| 5° p> Yo Mute 
m n 1 ,2) 1/2 
” 2 |All |S>13> FS Mt } , 








t=1 \j=1 


where ¥v,, and z, denote scalar components of v, and z, respectively, 
and A/J], denotes the i-th row vector of A/J. Hence for arbitrarily 
given positive numbers ¢ and 6 


Pr {| 5 vAllz, 





< 8, for |A/7| = Z 
(8) 


n 


= Pr[ ye 


t=1 


1 
m7 = 8 


Ni 








i” < | ; 


On the other hand, we have 


(8.1) E(U4 25 |2ns Zn, °°) = 0, 

(8.2) Var (4.2 |Z, 2, ***) << OF < @, 

(8.3) conditional random variables v,,z, given (Z,,2Z,.,°**) are in- 
dependent of each other with respect to different t, 
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from assumptions (ii), (iii), (iv). Therefore, we get 


N ) 
Pr {| Sot <%f BES Me (i = 1, +++, m; 9 =1, +++, n) 


for any positive 5,,»,, and any integer N=WN,, taking a sufficiently 
large integer N,, as a result of Bernoulli’s law of large numbers. 


Consequently, 
(9) Pr| 33435 | 2-2 
i=1 jJ=1 NG — 


for any positive 5, and \,, and N2WN,. Combining (8) with (9), we 
get the inequality 





2)1/2 
} <8,)21->, 


< 8, for |Al| =e} >1—2 








(10) Pr {|2-yo o Alla, 


4 


for any positive 6 and X, and N2WN,, taking a sufficiently large N,. 
As for last term in (7), we have following: 


7b Ale 


a at |z,|*| A/T, |? cos? (AZZ, 2,) 


z(- 


We (Allz,)' (Allz,) = 


Il 


a = = 


t=1 


a 
m,? + wlll cos* (A/T, 2) , 


where cos (A//,, z,) denotes the cosine of angle between vector A//, 
and z, in n-dimensiona! Euclidian space. Since the inequality 


lim +5 z,|° cos* (All, z,) 


im | P lim Ecos" (All,,z, =, >0 


N= 


= 


2 hi 


holds for A//J, such that |A//,| #0, from assumptions (v) and (vi), we 
have 


lim 5 > (Allz,)'(Allz,) = f(A) > 0 


as long as |A//| + 0, and, considering the continuity for f(AJ//) with 


respect to A//, 
inf f(All)=7>0 


[di|=e# 


for any positive ¢. Consequently we conclude that 
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(11) lim 5, 3 (Allz,)(Allz,) = 7 > 0 


for |A/7| = «. 


We now obtain the following result from relation (7), taking & in 
(10) smaller than 7 in (11), and combining (10) with (11): 


(12) Pr{SHYy, 1) > Si(Yy, 1,), for | — M,| =} 21—2r 
for any positive \ and «, and N2N,, taking a sufficiently large N,. 


If we substitute /] + (k — 1)AZ] or 17, + KAI] into Sj(Yy, 17) in (7) 
instead of /7, where k is a scalar, we get 


Si(Yy, 1, + kAll) 
= S3(Y¥q, I.) — 2 9 wikis Me, + 2. 5 (he 2,) (KIT) 


t=1 (=) ~ 


N N 
= Si(Yy, 1.) — = ky> viAllz, + & > (Allz,)'(Allz,) 


t=1 t=] 


= Si(Y¥q, Ml) + kL > walle, + <-Y yall \(All2,} 


t=1 


k(k —1)2 , 
+ es ame? (Allz,)'(Allz,) . 


Therefore, if we make k = 1 and are concerned with a point in Y, 
space such that the inequality 
Si(Yy, 1) > Si(Yy, My) 


or the equivalent inequality 


FE vAllz, + (Allz,)'(Allz,) > 0 


t=1 
holds, then inequality 
(13) Si(Yy, 1 + (k — 1A) > Si(Yy, 1) 


must hold for the same point of Y, space and the translated point 
of J], which can be written J] + (k — 1)(/] — I/,). However, any point 
in /7 space for which |/7, — /7| >« can be expressed as /], + (k — 1) 
(Ui, — 11,), using a point JJ, such that |/7, — /],| =¢, and a scalar k 
such that k=>1. So if the inequality 


Si(Yy, M1) > Si( Yy, I1,) 


holds for every point in a certain set of Y, space, and for every 
point of 77 space such that |/J — /J,| =«, then the same inequality 
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must hold for every point in the same set of Y, space, and for every 
point of /7 space such that |/J — /],|=>¢. Asa result we get 


Pr {Si(Yy, ll) > SiWYy, I1,), for | 7 = IT,| = &} 
= Pr {Si(Yy, 11) > Si(Yy, 11,), for | — 11,| >}. 

Thus we derive the conclusion from (12) and (14) that 
(15) Pr{S3(Yy, 11) > Si(¥y, 1), for | — M,| =e =>AQ—»), 
which holds for any given positive \,«, and N2WN,, taking a suffi- 
ciently large integer N,. 

In another context the inequality (15) can be transformed into 

Pr {| ly — 11,| <e} 21-2, 

because the fact that the inequality Si(Y,, 17) > Si(Yy, 17,) holds for 
all 17 satisfying |// — //,|=«, for a particular point of Y, space, 
implies that every point of /J space satisfying |/] — 17,| =¢ cannot 
qualify as //,, given the point of Y, space which occurs. 

From assumption (i), /7, and the consistent estimate //, of J],, ob- 
tained under the restrictions /J,¢€2,, can be uniquely transformed 
into B,, f°, and the estimates of B,, /, respectively, through a priori 
information about B and 7. Accordingly, estimates thus transformed 
necessarily become consistent estimates of the original structural pa- 
rameters, since this transformation must be continuous with respect to 
IIT, because of its uniqueness and the continuous mapping of B and 


I into 7. 

Finally, it may be remarked that if one is interested in the theorem 
only, the proof can be substantially simplified by starting from the 
consistency property of simple least squares estimation of the reduced 


forms. 


(14) 
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PROBLEMS IN PRODUCTION-INVESTMENT 
PLANNING OVER TIME* 


By VERNON L. SMITH’ 


1. INTRODUCTION 


IN A PREVIOUS PAPER [5] I attempted to develop the foundations for 
an integrated theory of production and investment in which the choice 
of best input-output technique, including the optimal stock levels of dura- 
ble inputs, and the optimal ‘‘replacement’’ of durable inputs were 
simultaneously determined in a cost minimization model. Very little 
was done in that paper to develop the dynamics of the theory, which 
is the objective of the present paper. 

We imagine individual decision making units, which we call ‘‘firms,’’ 
whose input-output decisions are constrained by a convex function re- 
lating the stock and flow inputs to a producing process. This function 
involves both flow and physical stock inputs because the physical stock 
of the typical durable input to production directly influences the trans- 
formation function of flow inputs into flow outputs.’ A most important 
characteristic of such individual durable goods is that they can be as- 
sumed to be continuously variable in size when considering expansion- 
ary investment decisions but cannot be contracted in amount, except 
by shutting down individual physical units. Typical examples are 
boilers, buildings, pumps, and pipe lines. For example, a steam gen- 
erating plant can install any size boiler, X, measured, say, in square 
feet of heating surface. But once a unit of size X° is installed, dis- 
investment can occur only by removing the unit entirely from produc- 
tion. Once a unit has been installed, expansion can take place either 
by adding a parallel unit of appropriate size or discarding the existing 
facility and replacing it with a larger unit. Hence, the decision prob- 
lem contains a curious mixture of continuity and lumpiness. 

The analysis will be confined to processes which require one current 
input and one capital equipment input. 

* Manuscript received January 15, 1960; revised April 26, 1960. 

' The research underlying this paper was supported by a Ford Foundation Faculty Re- 
search Fellowship, 1958-1959. I am indebted to the referee for his numerous detailed 
suggestions and corrections which have materially improved the mathematical exposition 
and analysis of the paper. 

* Ibid., |5, pp. 62-66}. 
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2. THE SINGLE INVESTMENT 


Suppose the firm’s planning horizon is 7 years. Assume that a 
single investment decision is to be made over the planning horizon, 
and that decision concerns the initial stock of capital to be installed 
by the firm. Assume further that output requirements as a continu- 
ous function of time are known with certainty, say y(t)0<t< T. 
The firm must choose the level of investment outlay to be made 
initially and the rate of consumption of current input at each instant 
over the horizon. The consumption of the current input at time ¢ is 
2,(t), while the quantity of the capital good to be installed is X,. The 
process is constrained by a production function of the form y(t) = 
f(z,(t), X,] which exhibits non-decreasing returns to scale. It is 
assumed that no output can be produced without some positive amount 
of each input. If the firm desires to minimize discounted total current 
account cost’ over the interval 0<t< T, then the problem is to 
choose 2,(t) > 0 and X, >0 so as to minimize 


MNT Y 3 -rt r WX, 
(2.1) C= |wn(tyerat + 7. , 
subject to 
(2.2) y(t) = f[x,(t), X;] , 


where w, and W, are the prices of the current and capital inputs, and 
r is the rate of interest. 
An equivalent problem is to minimize 


2.8) p=" ("waite — My latd), X) — yop yae + PMR 


2 
—e 
with respect to the extremal x,(t) and the variable X,, where X is a 
function of w,, W,, T, r, yi(t), and t, to be determined shortly. The 
function » plays a role entirely analogous to the Lagrange multiplier 
of the ordinary theory of relative optima. 

Note that the expression for ¢ in (2.3) is a special form of the 
following standard calculus of variations problem: Minimize 


8 “Discounted total current account cost,” or simply “current account cost,’’ as I call 
it, is the constant outlay stream that has a present value equal to that of all actual outlays 
for current and capital inputs over the horizon. For the models in this paper, minimiz- 
ing current account cost gives the same decision rules as minimizing the present value 
of all outlays. 
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(2.8’) Ie (Fe, é, thdt , 
where z = x(t) and 4 =dz/dt.* The solution to (2.3’) is contained in 


the Euler equation oF os (46 ) = 0 in x(t). We apply this condi- 
tion directly to (2.3) in minimizing ? with respect to the extremal 2,(¢), 
and we apply the condition 0~/0X, = 0 in minimizing with respect to 
the variable X,. If we let F’ stand for the integrand in (2.3), necessary 
conditions for an interior minimum of # are’ 


(2.4) OF — we- — rf{[x(t), X,] =0, 
On, 
PU Py =. / ey 
(2.5) Ox, = ee —| Vila, X,)dt + i_-ew ’ 


where f, = 0f/0x, and f, = 0f/@X,, and a similar conventional notation 
is used throughout this paper. 

From the production function (2.2) and the necessary conditions (2.4) 
and (2.5) we determine (if they exist) the optimal consumption path 
of current input x(t), the optimal level of investment, X}, and, as a 
by-product, the function \»° = (w,, W,, T, 7, yj(t), t) is determined. 
Specifically, by eliminating \ from (2.4) and (2.5), we get 


r w,f.{2,(t), X;] e~* rW, = 
e® | -o* \ f[a,(t), X,) th 1-o” f 


The interpretation of this expression is easily clarified. The ratio 
f./f, is the amount by which the consumption of input number 1 is 
decreased (increased) at time ¢ if the initial installation of capital is 
increased (decreased) infinitesimally. Then w,f,/f, is the rate of sav- 





* See, e.g., Fox [2], pp. 3-34. Also see Fox [2], chapter IV, passim, especially pp. 
80-83 and pp. 94-96, for a discussion of constrained extremal problems where the sub- 
sidiary equation is of nonintegral type. The problem presented by Fox (in one variable 


form) is to maximize or minimize = ["Gce, @,t)dt subject to S(z,z,t)=0. He then 


introduces 4, a function of t to be ae and develops a solution (ibid., [2], p. 96) 


of the Gian Se a) )- {a= 5S _ (2 
oz 
oF 


oz dt 
a 
tion oF (4) with F = G(z,2z,t) — AS(z,z,t). The well-known Lagrange 


multiplier shadow price tricks of static theory apply to our class of calculus of variations 
problems. 

® Note that the condition (2.4) can be obtained by straight scalar optimization, since 
the integral in (2.3) does not contain any time derivatives of the extremal. The calculus 
of variations form of the solution is emphasized because it underlies all the problems to 
be discussed, and some of the subsequent problems cannot be treated by scalar techniques. 


=} <0, This is equivalent to the solu- 
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ings in the cost of the current input per unit increase in the initial 
stock of capital. Hence, the integral of (w,f,/f,)e-" over the interval 
0<t< T is the present value of the savings in current input cost 
per unit increase in capital investment. Equation (2.6) states that in- 
vestment is expanded, given the time pattern of output requirements, 
until the discounted annual value of the current input cost saved from 
the last increment of capital equals the unit rent of that increment 
of capital, rW,/(1 — e-"7). Now, by integrating (2.6) in combination 
with (2.2), we determine the optimal values x?(t) and X3. With x(t) 
and X3 determined, (2.4) provides the Lagrangean function \°. Hence, 
an interior solution to (2.2), (2.4), and (2.5) can be written 


(2.7) xi(t) = h,(w,, W,, T, r, yi(t), t) , 
(2.8) X; = Aw, W,, yA yi(t)) , 
(2.9) = Mw, W,, T, 7, yi(t), t) . 


The similarity between this solution and the analogous static solution 
will be evident. 

Can we carry the analogy with static theory a step further and in- 
terpret \ as a marginal cost function? The answer is ‘‘yes’’. As is evi- 
dent from equation (2.4), when put in the form \ = w,e-"/f,[z,(t), X,], 
i is the discounted instantaneous (or short-run) marginal operating cost 
function. The effect of changes in w, or W, on the consumption of the 
current input x(t) and the level of investment X? parallel the effect 
of similar changes in the analogous static model. The higher the price 
of capital relative to the price of current input, the lower is investment 
and the higher the consumption of current input. 


8. THE TWO INVESTMENTS CHAIN 


Consider a slight extension of the problem discussed in the last 
section. Assume again a finite horizon of T years. with a known 
dynamic requirements function, y(t), defined over the interval 0<t< T. 
We divide the interval into two parts of unknown lengths, 7, and 7,, 
thus forming two operating periods, I and II, as shown in Figure 1. 
It is assumed that at most two investments are to be made in the 
interval 0<t< T. The first investment is at t=0, the second at 
t = T,, to be determined. We wish to determine conditions for the 
following: (1) how much capita! to install initially; (2) how much capital 
to install at 7,; (3) when to install (if at all) the second capital facil- 
ity, i.e., the determination of 7,; (4) the optimal instantaneous operat- 
ing rules in periods I and II; (5) the optimal consumption path for the 
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current input in periods I and II; and (6) rules for discarding either 
capital facility in period II if such action is optimal. To continue with 
the earlier example, we want to know how large a boiler to install 
now, and how large a boiler to install at any later time. We also 
want to know how to load the boiler or boilers at each instant over 
the interval, how much current input (fuel) will be consumed at each 
instant, and when, if at all, to discard either boiler. 

We assume momentarily that neither facility is to be idled tempo- 
rarily over any subinterval. If at any time it pays to withdraw one 
facility from production, it is assumed that the action is permanent 
and the facility is discarded at zero salvage value. 

Using a current account cost criterion, the objective is to minimize 


@1y ~~ pe { |i menerat + WX} 


+ {| ee + axolerat + Waka}, 


where X,, is the quantity of capital input installed at the beginning 
of operating period I (facility number 1), X,, is the amount installed 
at the beginning of period II (facility number 2), zi,(t) is the consump- 
tion path of current input by facility 1 in period I, 73,(t) is the con- 
sumption path of current input by facility 1 in period II, and 273,(t) is 


a aati lle 


I il 


Investment = Xz Investment =X2, + Xz 








qT T2 














FIGURE 1 
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the consumption path of current input by facility 2 in period II. The 
expression (3.1) is to be minimized subject to 


(3.2) y(t) = f'[xi.(t), Xa], 0O<st<T,, 
and 
(3.3) y(t) = f'[xi,(t), Xu] + f'[zi(0), Xn], r<st<T. 


The functions f' and f? refer to the production function for facilities 
1 and 2. That is, the production technology ruling at the beginning 
of period I[ is /f’[zx!,(t), X,,], while the production technology available 
at the beginning of period II is f*[x,(t), X.,.]. By time 7,, technological 
improvements in equipment and operating technique may have rendered 
f? entirely dissimilar to f'. We allow for this by permitting the pro- 
duction function to alter from one investment point to the next.° 

By introducing two Lagrangean functions \, and X., one for each of 
the constraints in (3.2) and (3.3), a problem equivalent to the one just 
stated is to minimize 


= {| mane — aise, Xul— ue)at + W.Xu} 
=e J0 


(3.4) — rl, (w,fai,(t) + ai(t)'e* — ul PLei(t), Xu) 
1—e" Vr, 


+ f*[x3,(t), Xn] — y(t)})dt + W.X,e} 


with respect to the nonnegative extremals 2},(t), x7,(t), z7,(t), and the 
nonnegative variables X,,, X;., T, (since T is fixed, T, is determined 
once we know 7;). 

Assuming an interior solution, for the present, the necessary in- 
tratemporal conditions for minimum #Y with respect to the current 
input consumption paths -',(t), zj,(t), and i,(t), in that order, are 


(3.5) we" — ru, fi[zi,(t), Xn] = 0, 0<t<T, 
(3.6) we" — AS i[zi.(0), X,,] = 0, 

T, t y T; 
(3.7) we" — rf i[zi,(t), X»] = a} stsiit 


The intertemporal conditions on the investment levels X,, and X,, are 


6 A more general way to introduce technological change is to write the second period 
production function in the form f%Azi2(t), Xx, T;), i.e., the ‘“‘state of the arts’’ is a func- 
tion of the point in calendar time at which the capital equipment of a process is installed. 
The class of problems under discussion is tedious enough without pursuing maximum 
generality 
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an — {lm slen, Xaldt 


+ \: rafilei(t), Xldt — W,} =0, 


T 

G9) = --_*{|"asiteh@, Xaldt — We} = 0. 
1 —e-* (Jz, 

Finally, the optimal timing of the second investment, T7;, must satisfy 

the condition 09/0T, = 0, or 


(3.10) w,xi,(T,) — w,[vi(T,) + «i(T,)] — rW,.X, =0. 


The conditions (3.5)-(3.10) contain practically the whole of a modern, 
marginalist, dynamic theory of production and investment. », in (38.5) 
and (3.8) is the discounted marginal operating cost function in period 
I, ie., \, = wie~"/f! over the interval 0<t< T,. X, is the discounted 
marginal operating cost function in period II. From (3.6) and (3.7) 
we have ’1,= we"/fi = we"/f?, for T,<t< T, during operating 
period II, i.e., the two facilities are operated at output (and current 
input) rates that equalize their instantaneous marginal costs (and mar- 
ginal productivities). 

In interpreting equation (3.8), note that the first integral under the 
brackets gives the present value of the period I savings in current 
input resulting from an infinitesimal increase in the initial investment 
level X,,. Similarly, the second integral gives the value of the corre- 
sponding period II savings. Therefore (3.8) requires the initial invest- 
ment to be expanded until the discounted annual value of the savings 
from the last unit of capital equals the unit rent of capital, r W,/(1—e-"”). 
Equation (3.9) requires the level of additional investment at 7, to be 
expanded until the annual discounted value of the savings resulting 
from the last unit of that capital is equal to the unit rent of capital. 

The first term in (3.10) is the instantaneous operating cost, at T,, 
of operating facility number 1 alone. The second term is the initial 
cost of operating facilities number 1 and 2 in parallel after the latter 
has been installed. Therefore, the difference between these terms is 
the initial net rate of operating cost savings resulting from adding the 
new facility. The last term in (3.10) is the instantaneous interest on 
the investment in a new facility. Therefore (3.10) says to time the 
installation of a second parallel facility when the rate of operating 
cost savings equals the instantaneous interest on the new investment. 
In a word, a new facility should be installed whenever the first ‘‘year’’ 
savings effected by the new unit will pay the first ‘‘year’’ of interest 
on the new investment. 
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So far the analysis has assumed an interior solution. Can we re- 
place the strict equalities in (3.5)-(3.10) with inequality signs and give 
Kuhn-Tucker interpretations to the surplus fulfillment of one or more 
of these conditions? We can (see, e.g., Morin [4], p. 448), and the re- 
sulting conditions will provide rules for an optimal equipment discard 
policy. 

Consider first the conditions (3.5)-(3.7), which we rewrite in the form 


(3.5’) we" eg MS ilzi,(t), Xx] = 0, 0 < t < ae 
(3.6’) wen" — ufi[zilt), Xu] = 0, 
(3.7) we" — rf ilz,(t), Xa] > at 1sts tf. 


If ‘‘>”’ holds in (3.5’) for some t = t,, 0 < t,, < T;, then 2},(t,,) = 0, i.e., 
at ¢,, it pays to withdraw the first facility from production. If this 
occurs in the first operating period it means that output requirements 
have declined to zero. If ‘‘>’’ holds in (3.6’) for some t =t,,, T, <ty< T, 
then 2i,(¢,.) = 0, and at ¢t,, the first facility is withdrawn from produc- 
tion in favor of the singular operation of the new facility installed 
at 7T,. Or, if ‘‘>”’’ holds in (3.7’) for some t,, 7, <t. < T, then 
xi,(t,,) = 0, and at t,, the second facility is withdrawn from production 
leaving the original facility only. Hence, by proper interpretation of 
the corner solutions implied by (3.5’)-(3.7’), we can obtain all possible 
equipment discard alternatives contained in the planning interval 0 < 
t<T. Note the strategic role played by the inequalities; they define 
implicitly the firm’s optimal equipment discard policies. Our formula- 
tion of the problem allowed explicitly for equipment purchases only, i.e., 
at 0 and 7;. The discard points, if they exist, are then evolved by inter- 
pretation of the equilibrium condition. Also note that the extension of 
our analysis to dynamic problems has not destroyed the intuitive useful- 
ness of the Lagrange multipliers; they continue to provide valuable 
service as shadow price evaluators. In this sense the conditions 
(3.5’)-(3.7’) require each facility to be operated at all times so that 
the internal instantaneous value of a unit of the current input, e.g., 
rs i[xi.(t), X,,Je* in (3.5’), does not exceed the external cost, w,, of 
acquiring a unit of this input. 

Similarly, we can replace ‘‘=”’ with ‘‘=>’’ in (3.8) and (3.9) to ob- 
tain more general conditions on the investment levels X,, and X,,. If 
**>”’ holds in (8.8), then X,,=0. This could only obtain if require- 
ments are zero throughout the first operating period, since we assume 
that no output can be produced without some positive investment. If 
‘*>” holds in (3.9), then X.,=0, and the best policy is to employ 
singly the original facility, X,, > 0, throughout the planning period. 
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The condition (3.10) could also be written with ‘‘>’’ in place of ‘‘=’’. 
If ‘‘>’’ holds, then 7, = 0, and only a single facility is used over the 
entire period. Hence, X,,=0 and 7, = 0 have equivalent implications. 

What happens if we drop the assumption that neither facility is to 
be idled temporarily over any subinterval? We can still write conditions 
on optimal policies, but the problem becomes more cumbersome. Now, 
if ‘‘>”’’ holds for some ¢ in any of the conditions (3.5’)-(3.7’), the cor- 
responding facility is discarded only if there is no later t’ > ¢ at which 
the ‘‘=’’ in that condition can be satisfied. For example, suppose 
“>” holds in (3.6’) for some interval T, < t,, <t< ti, < T, with the 
the equality holding elsewhere. Then at ¢,,, instead of discarding the 
equipment, it is placed in standby until ¢/,, at which time it is returned 
to the operating state. 


4. THE MULTIPLE INVESTMENT CHAIN 


Turning next to a generalization of the two-input problem in pro- 
duction-investment dynamics, our objective is to extend the analysis 
of the last section to the infinite horizon investment chain. We divide 
the planning interval, 0 < t < ~, into an infinite sequence of segments 
of unknown lengths 7,, 7,, 7;,++-, forming an infinite sequence of 
operating intervals, as shown in Figure 2. The objective in this case 
is to minimize the discounted current account cost over all future 
operating intervals, subject to input-output constraints of the form 
(3.2) and (3.3) for each operating interval. By introducing Lagrange 
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multipliers \,, \., 3, +++, this constrained minimum problem is equivalent 
to minimizing the tedious expression 


any &* re {| 2 (wx wiftje"* — hel OP TRO) z.j- vit)}) dt 


t=1 * j=1 


+ W,X0-7s-s} , 


where Tf = }0}., T,, and T} = 0, with respect to the extremals x(t), 
the capital stock levels X,,, and the purchasing dates T3(j < k =1, 
2,3,+++). The expression (4.1) is the infinite dimensional counterpart 
of (3.4). We again assume that whenever optimality requires a facili- 
ty to be idled, that facility is discarded. 

With respect to the extremals, necessary conditions for minimizing 
(4.1) are 


(4.2) we" — rf i[zi(t), X,] = 0, k=1,2,3, +++; 
41=1,2,-+-,k. 


The infinite dimensional conditions on the X,, and T, are 


(4.8) -r | 


rT, : 
[ir rusted, X,)dt} + reraW, > 0, 
bay ‘ST e-1 
j a 1, 2, 3, —— 


and 


oo & oo k 

p w.d0 in( Ti)e"" _ y Wd) inl Ti_e*e-1 
(4.4) kaj hel b k=J+1 heal 

—r >> W,Xye-7-1 > 0, j =1,2,38, ++ 
k=J+1 
where T? = }OtL, T, and T} =0. The intratemporal conditions (4.2) re- 
quire the current input to each facility in all operating periods to have 
a marginal inputed value which does not exceed the market price of 
the input. When the equalities hold, these conditions also imply that 
all facilities in any period must be operated at output levels which 
equate their marginal operating costs (and marginal physical produc- 
tivities). If, at any instant in any period, an equality contained in (4.2) 
cannot be maintained, the appropriate facility is at that moment idled 
and, therefore, by assumption, discarded.’ 
The intertemporal conditions (4.3) require investment in a new facili- 

ty at each purchase date to be expanded until the discounted current 
account cost of all future savings in operating expenses, effected by 


7 The problem of temporary idling can be handled by seeking subintervals over which 
the appropriate inequalities hold in the conditions (4.2). 
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the last increment of physical investment, is equal to the annual rent 
of the capital good. The integral in (4.3) serves the function of 
cumulating these discounted savings in each future operating period. 
These savings are then summed to obtain their total. If the rent of 
capital exceeds these savings for all feasible investment levels (the 
‘*>”’ prevails), then the investment is not undertaken. It is entirely 
possible that these conditions will signal an end to the investment 
chain. For example, if optimality requires the nth facility to be the 
last in the chain, then “>” will hold in (4.3) for X,, k=n+1, 
n+2,n+8,+++. The same situation will also be reflected in (4.4) if 
there is no solution for T,,,, i.e., no new facilities are acquired after 
the nth investment. In general, the conditions (4.4) determine the 
length of the operating intervals T7,, T,, T;,-+-, and therefore the 
purchase dates T'*, T}, T}, ---. These conditions require a new parallel 
facility to be installed whenever the discounted cost of operating the 
existing complex of parallel facilities is equal to or exceeds the dis- 
counted annual operating and capital cost of an additional facility. 
Each year we imagine the firm comparing the cost of maintaining the 
status quo (i.e., continuing to operate the existing equipment complex 
without benefit of an additional unit) with the cost of adding a new 
unit, proper account being taken of the savings effected by the latter. 
The status quo is maintained an additional year if the discounted cost 
of making the addition exceeds that of not making it. Abrupt tech- 
nological improvements, reflected by a shift in new equipment produc- 
tion functions, are, of course, favorable to the early purchase of a 
new unit. 


5. EXTENSIONS AND MORE GENERAL MODELS 


a. Discussion. The solutions of the last section are quite general 
for the two-input case. To obtain more specific results requires more 
specific information concerning output requirements, price parameters, 
and technological change. For example, in Figure 3, assuming an in- 
finite horizon, the solution might require a chain of five investment 
expansions to be made at A, B,C, D, and E. The third new facility 
might be added at C, even though output requirements have risen only 
slightly since the second facility was installed at B, because of techno- 
logical improvements occurring after the second facility was purchased. 
The investment chain might end at E, after five facilities have been 
purchased, if output requirements are due to decline to a constant level 
as shown, and no important technological improvements are in store. 
At F, ¢,, years after the purchase of the last facility, the declining 
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output level may precipitate the discard of the first facility (or any 
of the other five facilities depending upon the state of technology at 
the time each was purchased and the scale of the unit purchased). 
Can we generalize our two-input results to multiple inputs? The 
extension is easy if we have multiple current inputs, with a single 
durable input. The solution is identical with that of (4.2)-(4.4), except 
that we must add conditions of the form (4.2) for each current input. 
Such additional conditions simply require every current input to be 
consumed at rates which equalize their marginal costs, i.e., the ratio 
of the price of each input to the marginal productivity of that input 
must equal the ratio for every other current input in each operating 
period. The case of multiple durable inputs is not so simple. It is 
easy to see intuitively that conditions like those of (4.8) must surely 
apply to every type of durable input to the process. That is, at every 
purchase date for each type of equipment, investment is expanded 
until the savings on both current and capital account of the last in- 
crement of the capital good are equal toits annual rent. But to de- 
duce such conditions in a form similar to (4.3) is most formidable, 
since the moment we admit more than a single capital good we lose 
the one-to-one correspondence between the equipment purchase dates 
T*, Ti, Ty, +++, and the operating intervals 7,, T,, 7;, +++. That is, we 
must now distinguish purchase dates for each equipment type. In 
the general case of » current inputs and m capital inputs, we cannot 
even write an expression of the form (4.1) without knowing how the 
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purchase dates for each capital input are ordered with respect to each 
other. Furthermore, the problem cannot be side-stepped by attempt- 
ing to by-pass the production function and formulating the investment 
problem directly in terms of a current account cost function. There 
still remains the problem of determining such a cost function, and 
defining the optimal policies which it contains implicitly. The hypothe- 
sis of a cost function always presumes a wise class of optimal decisions 
to have been made. Also, the problem is not solved by formulating 
it in terms of ‘“‘capacity’”’ in the sense of output. This is a most 
nebulous concept, and leaves unanswered the question of the connection 
between ‘‘capacity’’ and the stocks of equipment in kind which are 
the direct objects of adjustment in investment decisions. 

b. The Solution with Maintenance Outlays. There are, no doubt, 
many different directions in which the dynamic models which we have 
developed can be extended. One important element which has not 
been incorporated into the dynamic analysis is the cost of capital 
equipment maintenance and servicing. It will be assumed that the 
rate of increase in maintenance costs (per unit increase in the 
consumption of the current input) is an increasing function of the 
cumulative consumption of current input. I have argued elsewhere 
[5], pp. 78-80 that equipment maintenance cost can be expected to 
vary with the equipment’s cumulative consumption of fuel or energy.’ 
The latter is a direct measure of the total exposure of a machine to 
wear and tear. If we let U, = U,(t) be the cumulative consumption 
of the current input up to time i, then U, = [zde)de. Our assump- 
tion is that dM/dU, = m(U,), where dM/dU, is the rate at which 
maintenance and servicing expenditures increase with increases in 
current input consumption, and m(U,) is an increasing function. Hence, 
the time rate of expenditure for maintenance and servicing is U,m(U,), 
where 2, = U,=dU,/dt. Both the current rate of equipment utiliza- 
tion and its cumulative utilization contribute to current maintance 
outlays. 

Suppose we illustrate the effect of maintenance cost considerations 
on dynamic solutions by means of the model of Section 2 above. Total 
current expenditure, including maintenance and servicing outlays, is 
now w,U, + U,m(U,). Making the change of variable x, = U, through- 


8 Also see [6] pp. 26-28 for an empirical study of truck-tractor maintenance costs, in 
which maintenance costs per mile are found to be an increasing function of cumulative 
mileage. In the present context, it would be argued that a more direct explanatory 
variable is cumulative fuel consumption rather than cumulative mileage. 
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out, the expression corresponding to (2.3) to be minimized can be 
written 


P= aoa. U, + U,m(U,)je-* — a4 f(0,,1X,) — y(t)} dt 


(5.1) pain! 
rW,X, 
+ i_<° 
The integral in (5.1) is of the form (2.3) and is minimized with re- 
- OF d/0F 2 

spect to the extremal U,(t) if —_-———(| —_} = 0. 

pec i(t) i aU, SAS 0. In this case the 
derivative with respect to U, will enter the solution, whereas in our 
previous analysis the derivative with respect to the time rate of 
change of the extremal was identically zero. In this model we get a 
true calculus of variation problem, whereas previously the analysis re- 
duced to ordinary scalar calculus minimization. 

Equation (5.1) is therefore mirimized with respect to U,(t) if the 
following differential equation is satisfied: 


(5.2) UM er — Lit, + m(U le" — FU, X)} = 0. 
au, dt 
Integrating (5.2) and solving for » gives 


[w, + m(U,)Je-"* — [o.gn-erat + C, 
f(U;, X,) 


where C, is an arbitrary constant of integration. Since 


(5.3) Mt) = 





TOM n>, _ (77 OM TOM 
[ore at = | Ue dt + \'Uohe dt, 


(5.3) can be written in the form 
_- 
\, Ogre rat + [w, + m(U,)Je* + K, 


fAU,, X,) 


If at t = T, the termination of the planning interval, we impose the 
condition (7) = [w, + m(U,)Je-"7/f,(U,, X,), then K,=0. The La- 
grangean X(t) is again interpreted as discounted marginal operating cost 
at time t. Equation (5.4) says that in equilibrium X(t) is the ratio of 
the discounted cost of employing an additional unit of input number 
1 at time ¢ to the marginal physical productivity of ‘input number 1 
at time t. The discounted cost of employing an additional unit of input 
number 1 at time ¢ is composed of three parts. The first part, 





(5.4) Mt) = 
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\" Uom/oU,e-"at, is the additional cost of maintaining the stock of 


capital over the rest of the planning interval from t to T as a result 
of consuming an additional unit of input number 1 at t. The second 
part, w,e-", is the direct cost of buying the additional input. The 
third part, m(U,)e~", is the additional cost of maintaining the capital 
in the current period, i.e. at t, when an additional unit of current 
input is consumed at ¢t. (Remember that the consumption of input 
number 1 exerts both a current and a cumulative influence on main- 
tenance outlays.) 

Finally, by setting 69/0X, = 0, we get the condition 


rT “ ; A 
(5.5) pal, —MOPUO, Xodt + 2M = 0. 


Equation (5.4) and (5.5), together with the constraint y(t) = f(U,, X,), 
provide necessary conditions for determining 23, X3, and ° over the 
interval 0 <t < T. 

We have written the solution with maintenance outlays for the case 
of a single investment and finite horizon. The solution is easily ex- 
tended to a chain of investments over a finite or infinite horizon. 
Where the Kuhn-Tucker-Euler conditions require a facility to be idled 
(or discarded), we get 2,(t)= U,t)=0, and total operating cost 
w,U, + U,m(U,) drops to zero for the idled facility. 

c. Inventory Policy and Production Planning; An Extension of the 
Modigliani-Hohn Model. An alternative extension of our fundamental 
models is to drop the implicit assumption of a nonstorable product. 
This introduces the necessity for simultaneous treatment of optimal 
inventory policy and optimal production-investment technique. We will 
follow a continuous version of the Modigliani-Hohn [3] model’ in in- 
troducing inventory considerations into our analysis. Modigliani and 
Hohn assume a given marginal production cost function which depends 
only on the rate of production at time t¢. In our analysis the short- 
run marginal cost function will be a by-product of choosing the best 
input-output technique over the planning interval. Certain interesting 
results follow from this construction of the problem. 

Following Modigliani and Hohn, we assume a finite planning horizon 
of T years, and take as our criterion the minimization of undiscounted 
cumulative cost over this period. Suppose the amount of product in- 
ventory held at time t is H(t). Thenif Y(t) is the unknown cumula- 


® See Morin [4], and Arrow, Kariin, and Scarf [1], pp. 61-69, for continuous formula- 
tions of the Modigliani-Hohn model. 
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tive production up to t, and S(t) is the given cumulative saies up to 
t, we have H(t) = H(0) + Y(t) — S(t). That is, inventory at ¢ is in- 
itial inventory plus the excess of cumulative production over cumula- 
tive sales at t. We assume one current and one capital equipment 
input to the process. We let k be the constant marginal cost of hold- 
ing a unit of inventory per unit of time. 

Cumulative cost over the horizon can now be written 


(5.6) C= [on + k[H(0) + Y(t) — S(@))dt + W.X. 

The problem is to choose 2,(t), X,, and Y(t) so as to minimize (5.6), 
subject to the production function constraint 

(5.7) f{z(t), X,] — Y(t) =0. 

Also, we have the important nonnegativity condition on inventory 
stocks, 

(5.8) A(t) = H0) + YOH-—-SHS0, 

and the boundary condition requiring no inventory at the end of the 
planning period, 


6.9) ¥(T) = | flz(t), XJdt = S(T) — HO). 
Necessary conditions for such a minimum are 
(5.10) w, — rf,[z,(t), X,] = 0, 
(6.11) [-ile, Xdt + W,=0, 
and 
(6.12) z—-®>d0. 
dt 


It is assumed that equalities hold in (5.10) and (5.11), since inequalities 
here are not of interest if some of both inputs are required to produce 
any positive output. .» is the short-run or operating marginal cost 
function, assumed given by Modigliani and Hohn. 

The new center of interest is the Euler condition (5.12). Integrating 
this condition from t to t + 6, where 0<t+0< T, gives 


{eae — Mt + 0) +t) >0, 


or 
(5.13) kd + X(t) > Mt + 4). 
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The left side of (5.13) is the marginal cost of producing a unit at 
time t and holding it in inventory until time ¢ + @.. The right side 
of (5.18) is the marginal cost of producing a unit at ¢+ 06. Hence, 
(5.13) states that the marginal cost of producing a unit at any given 
time must not exceed the cost of producing it at an earlier time and 
holding it until the given time. 

If the inequality in (5.13) holds over some interval ¢t, << t+06< t,, 
then in that interval it costs more to produce a unit at ¢ and hold it 
until any subsequent time, than to produce the unit at the later time. 
Therefore, at any t in this interval, we should make Y(t) as small as 
possible subject to the condition (5.8). If inventories are zero at the 
opening of the interval, i.e., H(t,) = H(0) + Y(t,) — S(t.) =0, then 
the smallest we can make Y(t) is Y(t) = S(t)— H(0), and current 
sales are satisfied out of current production in the interval. If in- 
ventories are positive at the opening of the interval, i.e., H(t,) = 
H(0) + Y(t,) — S(t,) > 0, then we make Y(t) = Y(t,), scheduling a 
zero level of current production, until such time as the opening in- 
ventory is exhausted. Thereafter, the previous solution applies, with 
current sales met out of current production. Of course, over any in- 
terval for which current production equals current sales requirements, 
i.e., Y(t) = S(t), the inventory part of the model disappears. The 
solutions for x°(t), X3, and »°(t) are then obtained by simultaneously 


wrreig 2 Poh yds  ceniier 
catistying |’ “wheat = W,, Sit) = fle(t), XJ, and Mt) =. 


If the equality holds in (5.12) over some interval, the solutions 
xi(t), X%, X(t), and Y%t) are obtained by simultaneously satisfying 
(5.7)-(5.12). However, in obtaining these solutions, the latter condi- 
tions are separable in a most interesting way. Suppose, for simplicity, 
we assume the equality in (5.12) to hold over the entire planning 
horizon. By integrating (5.12), and evaluating the integration constant 
at t=O and t = T, we get kt + (0) — M(t) = A(t — T) + A(T) — X(t) = 0, 
which defines marginal production cost as a (linear) function of time. 
Hence, insofar as the short-run marginal cost function depends upon 
time, this relationship is determined by the differential equation (5.12) 
independently of sales requirements, S(t), and the production function. 
This is in sharp contrast to the no-inventory case in which z/(t) and X} 
are jointly determined by the dynamic marginal-rate-of-substitution con- 


T . . 
dition | (w,f,/f,)dt= W, and the technological constraint Y(t) = S(t) = 
Sf (x,(t), X,]. Only after x°(t) and X$ are so determined do we obtain the 
equilibrium dynamic marginal cost function from \(t) = w,/f[x,(t), X,]. 
In the no-inventory case the dependence of marginal cost on time is 
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obtained only after we know the optimal path of current input con- 
sumption, 2,(t). But the moment we introduce inventories, fix sales 
requirements, and allow output to be a decision variable, the variation 
in marginal production costs with time is no longer due to the marginal 
productivity conditions and the production function constraint. With 
inventories, production planning is constrained by the requirement that 
the marginal cost of producing a unit at any time be equal to the 
marginal cost of producing a unit at any earlier time and holding it 
to the given time. The consumption of current input over time, and, 
therefore, the marginal cost function over time, must conform to this 
condition. The structure of optimal production policy is altered there- 
by in a manner which may not be prima facie obvious. 

Assuming the equality holds in (5.12), we can write the simulta- 
neous conditions on 2,(t), X,, Mt), Y(t), and 2,(7), as follows: 








Wy Da W, a 
64H M- N+ FE, x) Fee.” 
Tw, f.[x,(t), Xs] 9, — 
(6.15) \ af lat = W., 
(6.16) ¥(T) = \'#la(@), Jae = S(T) — HO), 
(5.17) Mt) = k(t — T) + Fiat) XT 


(5.18) Y(t) = [ fl=0, X,at . 


Equations (5.14)-(5.16) determine jointly the variables 2x((7') and X}, 
and the extremal x(t). Using the solutions X} and 2(T), we obtain 
Xt) from (5.17) as a function of ¢ and all the parameters of (5.14)- 
(5.16). Finally, using x{(t) and X:, we obtain Y°(t) from (5.18). The 
same procedure applies for any subinterval over which the equality 
holds in (5.12). 

In comparing these results with those of the Modigliani-Hohn model, 
the important thing to observe is the different role played by (5.12). 
In models of the Modigliani-Hohn type, the marginal cost function 
» = MY (t)] is a given function of current output independent of in- 
ventory and production planning considetations. Under such conditions 
(5.12) places a constraint on Y(t) only, whereas in our model, in which 
x is dependent upon choice of technique, the inventory condition (5.12) 


places a constraint on both X and Y(t). 
Purdue University, U.S.A. 
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SIR ROY F. HARROD ON INTERNAL 
AND EXTERNAL BALANCE* 


By MASAHIRO TATEMOTO' 


1. THE DISCUSSION of possible conflicts between the policy require- 
ments of full employment (without inflation) and those of equilibrium 
balance of payments is an important application of the Keynesian 
theory of employment. This problem has been one of the most im- 
portant issues in postwar international economics, and many writers 
have been concerned with it. Among them, Sir Roy F. Harrod, in 
the new edition of his book,’ devoted an ingenious chapter to ‘‘the 
theory of balance of payments by considering what should be done in 
order not only to get an equal balance of external payments, but also 
at the same time to maintain full employment at home.’’ In this 
short note we shall put his theory in a diagrammatic form and make 
the point more clearly. 

He defines the employment equilibrium or internal balance as Y‘= Y', 
where Y* stands for aggregate demand and Y* for the ‘‘supply po- 
tential’’ or ‘‘full employment income.’’ In the case where aggregate 
demand is less than the supply potential, realized income, Y, is equal 
to aggregate demand, i.e., Y‘ = Y < Y‘, and there is unemployment. 
If aggregate demand exceeds the supply potential, realized income 
equals the supply potential, i.e, Y* > Y = Y‘, where the inflationary 
gap‘ is measured by Y* — Y*. 

* Manuscript received October 16, 1959. 

1 The writer has profited greatly from many helpful suggestions and comments made 
by Professors Michio Morishima of Osaka University, Anthony M. Tang of Vanderbilt 
University, and Harry G. Johnson of the University of Chicago. 

3 It was the late Professor Ragnar Nurkse who first pointed out the importance of this 
problem. See R. Nurkse, ‘Domestic and International Equilibrium,”’ in S. E. Harris, ed., 
The New Economics, Keynes’ Influence on Theory and Public Policy (London and New 
York: D. Dobson, Ltd. and A. Knopf, Inc., 1947), pp. 264 ff. Professor J. E. Meade also 
discussed in detail the ‘‘conflicts between internal and external balance”’ in his The Balance 
of Payments (London and New York: Oxford University Press, 1951), Chapter X, pp. 114 ff. 

3 Harrod, R. F., International Economics, Fourth edition, Chapter VII, (Cambridge 
Economic Handbook, Nisbet and Cambridge University Press, 1957). 

* This excess of aggregate demand over the supply potential is called ‘‘the quantum 
of unfulfilled orders.” If prices rise so as to reduce demand to the supply potential, 
the inflationary gap will be cleared at once. The classical school is concerned with this 


case. But it is probable that prices will not rise so much and that the inflationary gap 
takes the form of unfilled orders or waiting list. This is the Harrod case. 
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External balance is defined as the equality of exports and imports, 
X = M, where X and M denote the exports and the imports respec- 
tively. In the event that X = M, there exists ‘‘external disequilibrium 
or imbalance.’’ ® 

If Y*= Y* and X = M, the economy is in equilibrium in both re- 
spects. But there may be an employment disequilibrium, Y‘* = Y’, 
while there is an equal external balance, X= M. Or there may be 
an unequal external balance, X = M, even though the internal economy 
is balanced, Y* = Y‘. Thus, Harrod sets out the following four com- 
binations of imbalance and proposes corrective remedies for them.° 














Case Situations Remedies 
I Y¢< ¥* and X> Mere yet | Retlate 
Il Y¢> Y* and X¥<M omen Senet) Disinflate 
m | 70< > ant <a, Coomera ot | Belen totes sarant 
V | re> roms zou, lite at, |... | Reaeee aie meee 





2. It may be of some interest to reformulate Harrod’s verbal dis- 
cussions in terms of the ordinary 45° line diagram of national income 
determination. Instead of Harrod’s concept of national income, defined 
as the total sum of the expenditure on home-produced consumer goods, 
home-produced investment goods, and exports (excluding import con- 
tents), we shall here use a more popular definition of national expendi- 
ture, i.e., 


where C, J, and G represent consumption, investment, and government 
expenditure respectively. Figure 1 is a slightly modified 45° line 
diagram of income determination. The level of national income is 
measured along the horizontal axis; exports, X, and imports, M, are 
measured downward along the vertical axis, and the aggregate demand, 
C+I+G+X—M, upward. Suppose the aggregate demand is rep- 
resented by the line aa’; then the income, Y, corresponding to the 
intersection A of -the line aa’ and the 45° line stands for realized in- 


® It is assumed that there are no gifts and no spontaneous movement of capital. Harrod 
distinguishes the imports of consumer goods from those of investment goods, but here 
we use the aggregate category, “‘the imports’’, for simplicity. 

* Harrod, op. cit., pp. 143-152. 
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come at which total demand equates total supply. The supply potential 
or full employment income is represented in our diagram by the income 
level Y’, and there is unemployment, which is represented by the 
distance YY*. Yet the income level Y* which warrants the external 
balance (given by the intersection of the lines X and MM) lies in our 
diagraiu to the left of Y, and there is an external deficit represented 
by the vertical distance EF. Thus, Figure 1 is the graphical repre- 
sentation of Harrod’s Case III, where unemployment and external 
deficit are jointly present. Suppose some expansionist measures (fer 
instance, fiscal policy) to increase domestic expenditures (I or G) are 
taken; then, in our diagram, aggregate demand will be shifted upward 
to bb’, and the full employment level Y‘ will be attained. This will, 
however, lead to the widening of the gap between imports and exports 
(viz., from EF to E’F’). On the contrary, deflationary measures will 
succeed in recovering the external balance by shifting the aggregate 
demand downward to cc’, but it will lead to an increase in unemploy- 
ment (viz., from YY* to Y*Y‘). This is Harrod’s ‘‘doubly misfor- 
tunate’”’ case in the sense that expansion aggravates the external 
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balance, while contractior disturbs the internal balance. Therefore, 
the increase in aggregate demand should not be effected by increas- 
ing domestic expenditures I or G, but by increasing X or reducing 
M. The latter wiil have the effect of shifting Y* to the right. The 
appropriate remedy will, according to Harrod, be a reduction (either 
selective or general) of domestic factor rewards or a devaluation of 
the currency. This will, subject to the exchange stability condition, 
lead to an increase in X and a decrease in M and cause the required 
shift to Y*. ; 

3. The above diagram can readily be modified to represent the 
other cases. If the aggregate demand is shifted upward to dd’ 
through an increase in domestic expenditures, we will have the in- 
flationary gap represented by the vertical distance DP. This inflation 
will cause an increase in imports and a decrease in exports by mak- 
ing home-made goods more expensive than foreign goods. The line M 
will shift downwards and the line X upwards, as represented by the 
dotted lines in Figure 1. These two effects will lead to an aggravated 
external imbalance. This is Harrod’s Case II, for which disinflation 
is proposed as an appropriate remedy. But what happens in the op- 
posite phase when aggregate demand is shifted downward to ee’ 
through a decrease in domestic expenditures? At the realized income 
level Y’ we have unemployment and external surplus. So far as the 
symptoms of imbalance go, the situation gives the appearance of 
Harrod’s Case I. But it is not. The reason for this can be seen by 
considering the applicability of the remedy proposed by Harrod for 
Case I. The proposed expansionist measures applied to this case will 
eventually lead to full employment, but with external deficit. Put- 
ting it in another way, the expansion will only convert this type of 
imbalance to another kind, which is characteristic to Case III; thus 
expansion cannot be used as a trump card to wipe out the imbalance. 
The above case is quite different from Case I; we shall call it Case I’. 

But how can we represent Case I in our diagram? This can be done 
as follows. Assume a situation where the relative positions of Y‘ 
and Y* are reversed, as represented in Figure 2. If the realized in- 
come Y lies to the left of Y’ and Y* in Figure 2, we shall have 
unemployment YY* and external surplus EF again. This is nothing 
but Case I, because expansion here, as expected by Harrod, will have 
the effect of increasing employment and reducing external surplus. 
It will eventually lead to full employment with external surplus E’F’, 
rather than deficit as in Case I’ above. Incidentally, Case IV can 
easily be represented in Figure 2 as the case where the aggregate 
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demand exceeds the full employment income OY"(=BY") by BC. 

According to the opinions of government economists here, an actual 
example of Case I’ is provided by the Japanese experience in the 
1957-58 recession, which was characterized by under-utilization of 
capital equipment’ and external surplus. It may be considered that 
Y* lay to the left of Y* in this Japanese situation. Expansion was 
out of the question because it came shortly after the government 
adopted deflationary measures to correct the big external deficit (which 
appeared in the first quarter of 1957, when the economy approached 
the full capacity ceiling). A similar situation is likely to occur in an 
economy which depends largely on imported raw materials. 

Cases I and I’ above are twin brothers but have quite different 
dispositions of their own. As the reader observes, this difference 
comes from the different relative positions of Y* and Y*, reflecting 
the structural basis of an economy. The policy of expansion should 
not be applied to Case I’, since it will in the process convert Case 
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FIGURE 2 


7? In a country like Japan where labor is abundant relative to capital, the supply 
potential Y* is determined by the full capacity of capital equipment, but not by full 
employment. The excess of the supply potential over realized income is measured by 
the under-utilization of capital equipment, but not by the unemployment of labor. 
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I’ to Case III. Appropriate remedies for Case I’ would then be a 
combination of the two remedies proposed for Case I and Case III. 
These are a policy of expansion until Y’ is reached and, beyond Y*, 
a policy of devaluation or factor reward reduction. The author 
believes that it is very important to recognize the practical implications 
of the relative positions of Y’ and Y*.° 


Osaka University, Japan 


8 Professor H. G. Johnson has kindly pointed out that Harrod distinguished between 
“the main line of cure’’ and ‘“‘the possible trimming’’ operation in an article preceding the 
new edition of International Economics. For example, referring to Case I, Harrod sug- 
gests: “If the internal stimulus has to be very great. .., the balance of trade may shift 
unfavorably, and if... it becomes negative..., it may in due course become needful to 
trim the foreign exchange rate downwards. If...an external surplus remains, an upward 
trimming may be appropriate.’”’ He also suggests the trimming operations for each of 
the other three cases. Thus, it is clear that Harrod was quite aware of the possiblity of 
Case I’, However, in his International Economics, there is no distinction between a 
major and a minor remedy. See Harrod, R. F., “Currency Appreciation as An Anti- 
Inflationary Device: Comment,”” Quarterly Journal of Economics, Vol. LXVI, No. 1, 1952, 
pp. 102-116. 
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TARIFFS AND TRADE IN A TWO-GOOD WORLD* 
By W. M. GoRMAN' 


1. INTRODUCTION 


THIS SHORT NOTE is a sequel to ‘‘The Effect of Tariffs on the Level 
and Terms of Trade’’ [1], some of the results of which are extended 
to the case in which there are more than two countries. 

The elasticity of demand for the exports of a single country is dis- 
cussed in Section 2. The argument relies heavily on the assumption 
that there are only two goods, so the conclusion that this elasticity 
is likely to be high need not be taken too seriously. Both argument 
and conclusions, however, are so simple that they seem worth present- 
ing; they may, after all, form the basis for more useful work.’ 

The general results proved for a two-country, two-good world in [1] 
are extended to one in which there are more than two countries in Sec- 
tion 2. Assuming that their elasi:cities of demand for import were 
j,k, their marginal propensities to import 9,70, and that they levied 
tariffs at the rate of 100 ¢ percent and 100 d percent ad valorem, 
the terms of trade were shown in [1] to be approx'mately 


i - as po 


while the volume of trade, appropriately defined, was shown to be 
approximately 


(+e, 


where m, 6, t are the appropriately defined average elasticity, propensity 
to import, and tariff level.’ The assumption that there are only two 
goods seems likely to be less harmful in this case. 

This seems an appropriate place for me to say that I wish to 
withdraw the suggestion, strongly urged in [1], that we may expect 
trade elasticities to be very high in practice. The argument for that 

* Manuscript received January 4, 1960. 

1 ] am grateful to Mr. Esra Bennathan, Dr. Frank Hahn, and D- John Wise for their 
comments. None of them should be blamed for my inadequacies. 

2 Besides, they are needed in Section 2. At least one reader found them very surpris- 


ing, so, perhaps, others may too. 
8 Units being chosen so that both the terms and volume of trade are unity in the 
absence of tariffs. 
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proposition depended critically on the assumption that the country in 
question produced a perfect substitute for the good it imported. Pro- 
fessor L. R. Klein has persuaded me that this is unreasonable.‘ 

A fellow economist who read this paper was under the impression 
that its conclusions were valid only if each country’s share of the 
market were fixed. The paper would not he very interesting if this 
were 30, for it would imply that the various offer curves on each side 
differed only in scale and therefore had the same elasticity. I differ- 
entiate first and then put the observed market share into my formulae, 
so that they can be applied to the actual world where market shares 
are changeable.’ 


2. THE ELASTICITY OF DEMAND FOR A SINGLE COUNTRY’S EXPORTS 


Consider a world consisting of N countries trading in two goods X 
and Y. A, exports z, units of X in return for y, units of Y («= 
1,2, -++,”), and B, exports y, units of Y in return for 2, of X (r= 
n+1,n+ 2,+++,N).° A,’s elasticity of demand for imports is j,; B,’s 
is k,. 

Assume that each country’s trade is balanced and that the world 
price of X in terms of Y is P.’ Then, 


(%.) y=Pe,; y,=P2,; and y= Pz, 
where 
(2) s=}o2,=S2; y= Lu=Ly 


stand for the volume of world trade in X and Y respectively. We 
can therefore define 


(3) M=al/z=yl/y5 4 =2/e=y,ly 


unambiguously as A,’s and B,’s share in the world market respectively. 

We shall now seek the elasticity (with fixed tariffs) of demand K, 
for A,’s exports. Todo so, we displace A,’s offer curve’ and see what 
happens to the net foreign demand 7, = }>2,— }orw2%,, for A,’s 
exports. 


‘ In combination with Sir Donald MacDougall’s fine study of British and American 
exports [2], which I have recently reread. 

° Of course the reader may find some of the other assumptions unrealistic. 

* Throughout this paper A and i will vary over the A’s, r over the B’s. 

7 Indeed we make all the usual assumptions of optimum tariff theory—except, of course, 
that the tariffs need not be optimal and that there may be more than two countries. 

8 Possibly by allowing A; to levy a new tariff, possibly otherwise. Note that we are 
estimating the foreign demand for Ai’s exports and must therefore displace its own offer 
curve. 
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By definition,’ 


(4) k, = —P92,/2,0P; 9, = Poy,/y6P, 
therefore, 

(5) Péx,/6P = —k,x, = —p,k,2 by (3), 
and, by (1), (8), and (4), 

(6) Péx,/0P = dy,/0P — x, = »»(7, — 1)2. 


(6) gives the change in the supply of A,’s exports when their price 
increases by dP. 

The net increase in the demand for A,’s exports is 

0z,/0P = >> ox,/0P — S> ax, /aP, 
avi 

therefore, 
(7) — Poz,/aP = [j + k —-1— 2,9, — 1), 
where j = }>X,j, and k= }> yk, are the average import elasticities 
of demand for Y and X respectively.” 

Initially, Z, = x, = 4.2, so that 


(8) K,=G+k-—D/y,-— &, — D. 
The elasticity of demand for B,’s exports is similarly 
(9) J,=(9 + k—1)/p, — (&, — 1). 


Unless A, dominates the world market, the leading terms in (8) and 
(9) will be the more important, so that” 


(10) K,=(G+k—-—D/y and J, =(§@+k—1)/z, 
and, in particular, 
K,=J, if \=#. 


Hence, the elasticity of demand for a country’s exports depends 
mainly on its share of the market and is largely independent of the 
good which it exports. It seems likely to be high, too. If j,=j= 

9 i, = —d log y%/d iog (1/P) = Glog yi/Alog P. The partial @’s are used to denote the 
fact that tariffs are held constant in the other countries. The demand and supply re- 
lations are given in (11) below. 

0 Appropriately weighted. 

If all the j’s and k’s are greater than one: 


G-1- AG —- Y= Alfa - NSO, 
hi 


so that Ki 2 k/Ay and Jy = j/pr- 
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k = 2, for instance, and A, has twenty percent of the world market, 
K, = 14." 
What happens if other countries change their tariffs as well? 
Suppose A, levies a tariff of 100 t, percent ad valorem, (i=1, 2,-++,n), 
and B, one of 100 t, percent, ("= "+ 1,--+,N). We can then write 


(11) 2, = F,(P, t,), x, = GP, t,) 
and find 

(12) Pdz, = (j, — 1)adP + Pxidt,, 
(13) Pdz, = —k,x,dP + Pxidt,, 
where 

(14) x; = 6z,/dt,, x}, = 02,/dt,, 


with the aid of (4). 
Hence, by (2) and (7), 


(15) (j + k —1)adP = P(S> zidt, — > aidt,). 


Substituting into (12) we find the total elasticity of demand for A,’s 
exports, 


(16) K, = [649 +k -—DAJ—-G—- DV, 
where 
(17) 6, = xidt,/(S_ xidt, — > xidt,). 


If A, is the only country to change its tariff, 0,=1, and (16) is 
the same as (8)—as, of course, it should be. 

If A’s competitors—the other A’s—raise their tariffs to the same ex- 
tent, but the B’s do nothing, and if the A’s are similar to each other 
in the sense that 2,/z, is the same for all of them,” 0,=,, and 
K,=j-—i,+k. Finally, if all this is true and j,=j, so that A, is 
truly an ‘‘average’’ country"—except, perhaps in size: 


(18) K,=k. 


In other words, if the A’s are sufficiently similar” and act in a suf- 
ficiently similar manner, each of them is effectively in a monopoly 
position, as one would expect. 


In general, however, we know nothing at all about K,, not even 


13 Please remember that these results depend on the assumption that the consumer is 
indifferent between units of, for example, X made in different countries. 

8 Quite a reasonable assumption, since it refers only to proportional shifts. 

% Alternatively, one might argue that (18) were true ‘‘on the average’’ in some sense. 

% They may, of course, be different sizes. 
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its sign. The implications for the theory of optimum tariffs are ob- 
vious. On the other hand, we must remember that these results have 
been proved under the assumption that the goods exported by any 
one of the A’s are perfect substitutes for those exported by the other. 


8. TARIFFS AND TRADE 


In this section we shall show that the results proved in Sections II 
and III of [1] carry over to a many country world, at least as first 
order approximations. 

It will be convenient to consider first the effects of tariffs on the 
terms of trade P,” and to normalize our data so that, under free 
trade, 


(1) e=ny=P=1, 
(2.15)” implies that 
(2) dlog P = (3 xidt, — D aidt,)(j + k — 1). 


To convert this into a useful form, we must express xj and 2} in 
terms of the marginal propensities p, and oc, to import into A, and B, 
respectively. In order to do so, we first define 


(3) Q, = Pi(1 + t,) 
as the domestic price of X in terms of Y in A,, and 





(4) Z,=Y; — V0, = Q,2,t,, since y, = Px,, 
as A,’s exogenous income. (2.12) and (3) imply that 
(5) dy, = Pdxz, + «dP = (j,Q,”, + xidt, 
when Q, is held constant. At the free trade equilibrium, 
(6) Ps = (0Y9Z,)o, constams = Jy + Vil%. 
Similarly, 
(7) Oo, = k, + x/a,.” 
Substituting into (2), we find tha’ 
(8) log P = Ga pe  E—oK 


to the first order for ‘‘small’’ tariffs, where 


8 Strictly speaking, P is the terms of trade of the B’s. It corresponds to 1/P in [1). 
We refer thus to equations in other sections. 
8 This can easily be verified. 
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$=D xis k= Des; 





(9) p=»; =D pe,, 
e= Duis — Prt. d= L Uk, — o,)t, 
y i Mui, pane P) a LAK, ~— 0,) 


are appropriately weighted averages in each case. This is the result 
which we found in [1] for trade between two countries A and B, 
whose marginal propensities to import were @ and g, whose elastici- 
ties of demand for imports were j and k, and which levied tariffs of 
100c and 100d percent ad valorem respectively.” It would therefore 
seem that our two groups of countries can be aggregated to form two 
imaginary superstates whose elasticities, marginal propensities to im- 
port, and tariffs are ‘‘averages’’ of those for their component econo- 
mies. This sample pattern of aggregation is possible only because 
there is no trade between the members of the same group, and be- 
cause the tariffs are sufficiently small to ensure that no exporter of 
X becomes an exporter of Y, or vice versa. 

We can deduce from (8), as we did in [1], that tariffs are likely to 
have a negligible effect on the terms of trade.” 

Let us next consider the volume of trade, which we define as being 
I=V zy. Under our normalization, J=1 under free trade. 

Eliminating x and P from (2.3), (2.12) and (2.13), and using (8) and 
(9), we find that 


dlog I= dlogz + 0.5d log P, 
therefore, 
(10) log I= —(m — @)t, 
where m, 0, and ¢ stand respectively for the ‘‘average’”’ elasticity of 
demand in the export market, the ‘‘average’’ marginal propensity to 
import, and the ‘‘average”’ tariff, as defined in [1]. The averaging 
process is here carried out over all N countries. 

This again corresponds to the main result in Section II of [1]. The 
precise definitions of m,@, and ¢t are given in Appendix B to that 
paper. The quality of the approximation is only slightly impaired if 
we take 

m — 0.5 =the harmonic mean of j — 0.5 and k — 0.5, 
(11) @ = the arithmetic mean of ¢ and a, 
1+ ¢=the geometric mean of 1+ c and 1 +d. 


1° Of course, this is merely a first order approximation here. The result in [1j looks 
rather different, but is the same to the first order. 
* See the example in [1]. 
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The main theoretical results proved for a two-country, two-good 
world in [1] are seen to hold for a many-country, two-good world, 
even if in a somewhat attenuated form.” The main effect of the ex- 
istence of more than two goods is, presumably, the possibility of 
multilateral trade,” which is probably more vulnerable to tariffs than 
bilateral, because it has to surmount several barriers. These formulae 
may well understate the effect of tariffs in a many-good world, but 
it is not clear by how much. There is some discussion of this point 
in [1). 


University of Birmingham, England 
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% Though in an attenuated form: (i) these are approximations for small tariffs, (ii) 
there are no analogues of Propositions I and II of [1]. 

= Of essentially multilateral trade, that is, in which A sends X to B, B sends Y to 
C, and C sends Z to A, for instance. 








INTERNATIONAL ECONOMIC REVIEW 
Vol. 1, No. 3, September, 1960 
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1. INTRODUCTION 


RESEARCH IN long-run economic development is today one of the most 
important objectives in economics. Among the foundations for these 
investigations are the works of J. M. Keynes, who made the deter- 
mination and fluctuation of national income his primary goal. Since 
then, starting from the determination of national income and employ- 
ment in a closed system, models have been devised which make pos- 
sible an analysis of the process of economic growth. These models 
analyze many relevant variables such as income, employment, iavest- 
ment, and consumption as they change over time. 

With his theory of national income Keynes demonstrated that a 
whole economy may be represented by certain macro-economic aggre- 
gates, and that general economic equilibrium can be described by a 
few strategic variables. In this way, rather independent statistics of 
income, savings, and investment may be used to verify statistically 
economic theory. Keynes’ theoretical starting point, however, is the 
concentration of an infinite number of elementary economic units into 
“‘aggregates relevant for the functioning of the economy.’’' His ob- 
jective is thus a short-run analysis of the economic process ard not 
an examination of long-term growth. Keynes recognized these limita- 
tions and mentioned them explicitly in many places.? Similarly, many 
extensions of his work by his successors are limited to typical short- 
run problems. 

Many attempts have recently been made to set up theoretical 
models which will study long-run aspects of growth, utilizing the in- 
dispensable instruments of Keynesian theory. The analysis of short- 
run fluctuations around the trend of ‘‘equilibrium’’ is replaced by 
study of the trend itself. There had been attempts by classical econ- 
omists to understand long-run economic development, even if only in 

* Manuscript received February 25, 1960. 

1 Erich Schneider, Einfiihrung in die Wirtschaftstheorie, Vol. | (Tiibingen: 3rd ed., 
1955), p. 50. 

2 Gottfried Bombach, ‘‘Zur Theorie des wirtschaftlichen Wachstums,’’ Weltwirtscha/t- 
liches Archiv, Bd. 70 (1953, I), p. 110. 
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terms of secular stagnation; now, however, there is a renascence from 
the type of price and market theory—microeconomic, polypolistic, 
without regard to time—which for so long dominated economic think- 
ing. Harrod, Domar, Baumol, Hicks, and Bombach have devised im- 
portant long-run models; in particular, Harrod and Hicks have worked 
on a theoretical synthesis of short- ind long-run analysis. 

Although current investigations of economic growth are still predom- 
inantly theoretical, it has proved possible in some cases to use for 
verification statistical material originally prepared for the typical 
Keynesian short-run analysis. The study of changes in the level of 
national income over time has long been of primary importance in 
economic research, especially in growth analysis; empirical testing of 
models of economic growth should therefore start with the statistical 
determination of the nominal value of national income over a continu- 
ous long period. 


2. ORIGINS OF NATIONAL INCOME STATISTICS IN GERMANY 


In the case of Germany, quite extensive statistical material is availa- 
ble over a long period which, if properly prepared, could throw light 
on national income. This situation appears to be rather unique, since 
in other highly industrialized economies the relevant statistics do not 
go as far back into the last century as do the German ones. Now 
at last the scientific and methodological objectives of the German 
Historical School (which originally dominated German economic thought 
and later was almost entirely superceded by the Mathematical School) 
may be justified, ironically enough, by the very theoreticians who 
earlier fought it. The collection and publication of very detailed data 
about many different time series such as population changes, income 
tax payments, etc. are the results toa very great extent of the em- 
phasis which the Historical School placed on inductive methods. The 
importance of these series should not be overlooked, even though they 
were often too extensive and badly arranged and, in any case, be- 
came of real significance only after evaluation by deductive methods. 
The representatives of the German Historical School—like true children 
of their time—were to a large degree simply following a trend current 
in all sciences, natural, social, or humanistic, by attempting to push 
forward the boundaries of knowledge through the encyclopedic collec- 
tion of facts. 


8. THE HOFFMANN-MULLER STUDY OF GERMAN NATIONAL INCOME 


For these reasons, it has been possible to utilize the German figures 
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in order to verify certain aspects of modern growth theory. Professor 
Walther G. Hoffmann of the University of Miinster and Professor J. 
Heinz Miiller of the University of Freiburg, together with their staff 
of assistants and research workers, have succeeded in constructing 
for the first time national income figures for Germany over a period 
of more than one hundred years, i.e., from 1851 onwards. Moreover, 
these series are presented in a form which permits comparisons.’ 
These statistical results may make possible a threefold study of the 
process of growth of the German economy: first, the inductive explana- 
tion of growth; second, the testing of dynamic theories; and third, 
evaluation of future opportunities for development. They may also 
provide a basis for the determination of structural changes in the 
German economy from 1851 to 1957 by using cross-sectional stocktaking. 

Understandably, a collection of statistical series going this far back 
into the nineteenth century presents extraordinary problems and dif- 
ficulties. The authors must, however, be commended for taking this 
risk, for the results that have been achieved fully justify the under- 
taking. The book is limited to a reproduction of statistical material 
in the past hundred years; any theoretical interpretation of the re- 
sults must be left to further investigators, since a structural analysis 
of the sources of national income is not given.‘ The greater part of the 
Hoffmann-Miiller volume is devoted to an extensive examination of the 
nominal value of national income. Estimation of a parallel time series 
of real national income is hindered by the lack of price indexes proper 
for this type of macroeconomic investigation. Any price indexes 
available so far back usually omit one or another important component, 
the inclusion of which is essential for proper deflation of the nominal 
values. 


Let us take, for example, price indexes before 1918. Price changes 
in essential foodstuffs, which frequently move contrary to the cost of 
living index, are overrepresented; in addition, the available index of 
wholesale prices does not sufficiently reflect the development of pur- 
chasing power. Taking these difficulties into account, the authors 


5’ W. G. Hoffmann, J. H. Miiller, with the collaboration of Heinz Kénig and Klaus 
Déring, Jochen Dohmen, Rolf Goldschmidt, Franz Grumbach, Alexander Hélling, Arnold 
Hilden, Das deutsche Volkseinkommen, 1851-1957 (Tiibingen: J. C. B. Mohr (Paul Siebeck), 
1959). 

* Estimates for the generation of national product and the spending of national income 
are still lacking. The first problem will be extremely difficult to solve, but the latter 
has been attempted by Professor Hoffmann and Dr. Franz Grumbach ina paper presented 
to the Sixth European Conference of the International Association for Research in Income 
and Wealth, held at Portoroz, Yugoslavia, in August 1959 
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constructed from existing wholesale prices a new index which is made 
up of the individual indexes of wholesale prices for small groups of 
commodities, each weighted by its proportion to total household ex- 
penditure. Beginning in 1913, the official cost of living index is used. 
In this way the true growth of the economy has been indicated by 
deflating the money value of national income throughout the period 
under investigation. 


4, ESTIMATES OF GERMAN NATIONAL INCOME, 1851-1957 


The only areas in which an attempt has been made to determine 
national product generated have been agriculture and industry; the 
other sectors of the economy have been disregarded.’ Previously, 
German national income had been estimated on the basis of distribu- 
tion for the period from 1891 to 1913, although the Reich Statistical 
Office based its official study only on tax statistics of two member 
states of the German Reich.’ The same office computed annual figures 
for distribution of national income during the interwar period and 
increasingly comprehensive figures for generation of national product 
in agriculture and industry. Since 1948 the West German Federal 
Statistical Office has made annual estimates of national income and 
its most important components in the German Federal Republic. In 
1957 these figures were recomputed and the results used in the Hoff- 
mann-Miiller study under discussion.’ 

The methodological bases of the investigation are the concepts and 
definitions of income streams used in the standardized system of social 
accounts of the O.E.E.C.° In this system the principal streams of 
distribution of national income are: (1) employee’s compensation; (2) 
entrepreneurial and property income accruing to households and private 
nonprofit organizations; (3) corporate savings; (4) direct corporate 
taxes; (5) entrepreneurial and property income of government, less 
(6) interest on the public debt. The main currents of national income 
are estimated for the changing area of Germany during the entire 
period 1851-1957. In addition, the principal streams and their com- 


6 Cf. Paul Jostock, The Long-Term Growth of National Income in Germany, Series 
V, Income and Wealth (London: 1955), pp. 79 ff. 

® Das deutsche Folkseinkommen vor und nach dem Kriege, Einzelschrift zur Statistik 
des Deutschen Reiches, No. 24 (Berlin, 1932). 

7 H. Bartels, K. H. Raabe, O. Schérry, et al., ‘“‘Die Neuberechnung des Sozialprodukts 
fiir die Bundesrepublik Deutschland,”’ Wirtschaft und Statistik, N. F., Vol. 9 (1957), pp. 
123 ff. 

8 Organization for European Economic Cooperauon, A Standardized System of National 
Accounts (Paris: 1952), pp. 19, 75 ff. 
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ponents are given in detail for the following individual German states: 
Prussia, Saxony, Hesse, Baden, Wiirttemberg, Bavaria, and for the 
city-states of Hamburg and Bremen for the period 1851-1913. 

The statistical material available for the different streams of na- 
tional income is extremely heterogeneous. This is particularly true of 
income tax data, since until 1920 income tax legislation was left in 
the hands of the individual states. For this period, the computation 
of income streams for the entire German Reich is based on estimates 
for the individual states. For the period until 1870, indeed, only re- 
sults of Prussian income statistics are available; after that date, 
figures for other states are found increasingly. National income of 
the German Reich from 1851 to 1913 is estimated on the basis of per 
capita figures which are available for particular regions. It is assumed, 
obviously, that the results for the individual states are representative 
of the entire area; differences in income level resulting from the inclusion 
of data on federal states other than Prussia are adjusted after 1871. 
After 1920, national income is estimated on the basis of tax assess- 
ment data of the Reich. 

Tax collection statistics of the federal states until 1913, and, since 
1925, those of the Reich, are used to estimate entrepreneurial and 
wage and salary income of private households.’ Before 1891, when 
Prussia reformed its Income Tax Law, some components of income 
had not been included and made subject to tax; estimates for the 
period before 1891 had to be adjusted by an ‘‘Unterbewertungszuschlag,”’ 
i.e. an increase to take account of original undervaluaticn. Undistrib- 
uted corporate profits are not included in tax statistics and had to 
be estimated in special studies. In contrast to the system of stand- 
ardized accounts of the O.E.E.C., direct taxes paid by corporations 
are not separately computed to avoid double counting, since these taxes 
are paid out of undistributed profits and thus are included in the 
total figure. Data on government entrepreneurial and profit income, 
as well as interest on the public debt, are obtained from fiscal statistics 
and annual reports of social insurance agencies. 

During the past century Germany has undergone many territorial 
changes, so that it is not possible to use the same geographical area 
over the entire period under investigation. This, of course, makes it 
difficult to use different periods for purposes of comparison. 


® The period between 1914 and 1924 has been omitted; the wartime and post-war infla- 
tion do not yield reasonable figures which fit into the long-term trend. 
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5. RESULTS OBTAINED IN THE HOFFMANN-MULLER STUDY 


All these problems, which are here mentioned only briefly, are the 
subject of rigorous scientific criticism in the Hoffmann-Miiller study. 
This alone would render the book important, but it is, in addition, an 
excellent analysis of long-run national income computation. We shall 
not attempt to present all of the many rewarding results of the book 
but shall merely set forth a few of its important findings. 

The development of nominal and real German national income be- 
tween 1851 and 1957 is given in five-year averages, which affords a 
sufficient degree of precision for long-run analysis. The conclusions 
for single years may be dubious, but the correctness of the emerging 
trend is not affected. For an exact analysis of business fluctuations 
and cycles one must go back to annual values. These, however, are 
by their very nature not precise, since the most important components 
of income are considered in three-year-moving-averages in order to 
dampen oscillations and postpone cyclical turning points. 

The expansion of population throughout the entire period is covered 
by computing per capita income. From 1851 to 1957 nominal national 
per capita income increased by more than ten times the initial value. 
In order to estimate changes in real per capita income the authors 
use their newly constructed cost of living index to deflate nominal 
income. Even with this adjustment real per capita national income 
rose to four times the original figure. 

The period from 1851 to 1913 is of special interest, since it indicates 
the rapidity of economic growth with the increase in industrialization. 
In Germany the process of industrialization did not begin before the 
middle of the nineteenth century, while in England and France it 
started considerably earlier. The analysis in this book, therefore, 
provides a picture of industrial development almost from its very be- 
ginning. This is documented by the very low per capita income in 
the 1850’s. From that time on, the per capita income rose quickly 
and steadily. Between 1851 and 1913 the population increased by some 
80 percent, while the addition to total national income amounted to 
about 390 percent. That is, while population doubled and national 
income increased fivefold, the increase of per capita nominal income 
was about 170 percent, or almost triple. If one takes further account 
of the development of prices during this period, it is found that per 
capita real income is continuously above nominal income. Only the 
“‘Griinderjahre’’ (years of wild speculation) of 1871 to 1875 present an 
exception to this general conclusion. The period from 1851 to 1913 
represented a steady increase of real per capita income, partly because 
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of a rise in nominal income and partly because of falling prices (up 
to 1900). After prices began to rise, the purchasing power of nominal 
income decreased, but real income continued to rise. 

When we evaluate these findings, it must be remembered that 
changes in the relationship between the total population and the 
active working population have not been taken into account. Income 
per capita for the total population is, therefore, never identical with 
income per capita for the working population, nor is there a constant 
ratio between them, in view of the substantial increase in the ratio 
of the labor force to the total population since the first World War. 
The only conclusion that can be drawn is that real per capita income 
of the total population has roughly doubled in each sub-period of fifty 
years. 














TABLE 1 
NOMINAL AND REAL NATIONAL INCOMES IN GERMANY 
Nominal National Income Real National Income 
se Population , 
FRsuat | Grtiow | Taal, Yer copia [Race tndes| Taal, | Pee 
of marks, RM/DM) of marks, (marks, 
RM/DM) RM/DM) | RM/DM) 
1851-55 35,944 9,569 266 90.2 10,609 295 
1856-60 36,884 10,755 292 96.3 11,168 303 
1861-65 38,773 11,894 307 94.3 12,613 326 
1866-~-70* 40,141 13,051 325 96.7 13,372 333 
1871-75 41,641 15,171 364 103.4 14,672 352 
1876-80 44,105 16,280 369 92.1 17 ,676 401 
1881-85 46 ,039 17 ,557 381 82.0 21,411 465 
1886-90 48,176 20,104 417 81.3 24,728 513 
1891-95 50,825 22,638 445 80.2 28,227 555 
1896-1900 54,405 27 ,028 497 79.4 34,040 626 
1901-05 58,612 31,548 538 89.2 35,368 603 
1906-10 62,863 39,919 635 92.5 43,156 686 
1911-13 66,161 47,374 716 98.3 48,193 728 
1925-29 63,221 64,885 1,026 147.5 43,990 696 
1930-34 64,928 49,849 768 128.8 38,703 596 
1935-39 67 ,979 70,215 1,033 124.9 56,217 827 
1950-54 47 ,825 99 , 297 2,076 209 .0 47,511 993 
1955-57 49,825 148 ,277 2,976 221.1 67 ,036 1,346 


























Source. W.G. Hoffmann, J. H. Miiller, et al., Das deutsche Volkseinkommen, 1851- 
1957 (Tiibingen, 1959), p. 14. The territory covered changes. 
* Average of the years 1866 and 1869. 
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Cross section analysis reveals that private households’ share of total 
national income always amounts to more than nine-tenths, while 
undistributed corporate profits and government entrepreneurial and 
property income were quite insignificant at the beginning of the period. 
The percentage of national income going to the latter two recipients 
rose considerably, since their incomes grew more rapidly than those 
of private households. This is also true of interest on public debt, 
which in the process of industrialization increased extraordinarily with 
the greater role played by the state. This growth of interest on 
public debt was twice interrupted as a consequence of currency re- 
form at the end of each World War. 

Nominal income values are given by state from 1851 to 1913, but 
corresponding real figures cannot be computed because of lack of re- 
gional price indexes. In addition, the states investigated do not form 
definite and distinct economic regions, nor can they be characterized 
as predominantly industrial or agrarian. Very often the differences 
within a single state are of more importance than structural varia- 
tions between states. A more precise study of regional differences in 
welfare and wealth (real income) wi!l become possible only when small- 
er regional units are taken as the base for the analysis of national 
income, i.e., if the economic structure of the region is examined to 
detect differences in the composition of regional income. The very 
valuable study of Professors Hoffmann and Miiller and their collabora- 
tors is a first step in this direction. 
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SPRINGER-VERLAG IN WIEN 


Nationalokonomie 


Grundlagen fur eine exakte Theorie 
Von 
Dipl.-Ing. F. Korenjak 
Wien 


Mit 84 Textabbildungen. XXIV, 374 Seiten. Gr.-8°, 1960 
Ganzleinen $11.45 


Die theoretische Nationalékonomie ist noch weit vom Status einer exakten Wissen- 
schsft entfernt. Die vorliegende Arbeit hat sich die Aufgabe gestellt, ein Funda- 
ment fir eine cxakte Nationaléko-omie zu errichten. 
Die in fritheren Zeiten verschiedentlich aufgestellte eee daB die National- 
dkonomie yor eine exakte Wissenschaft werden kénne, wird wohl kaum mehr ernst 
ee es fehit nicht an Versuchen, cine derartige Theorie von verschie- 
Seiten in Angriff zu nehmen: Die »,Okonometrie‘* wht in analoger Rich- 
tung wie die exakten Naturwissenschaften vor, indem:sie versucht, von den Er- 
scheinun induktiv auf die Gesetze des Wirtschaftens zu schlieBen (wie der Autor 
einlei ausfiihrt, kann dieser Methode nicht der gleiche Erfolg beschieden sein 
wie in den Naturwissenschaften, da in der Wirtschaft keine Versuche mit Varia- 
tion von nur wenigen GréGen angestelit werden kénnen); rein mathematische 
stellen Axiomensysteme auf — eine Arbeit, die erst wertvoll wird, so- 
bald man wei, welche Axiomensysteme der theoretischen Nationalékonomie adiquat 
sind — die Wiener Schule schlieBlich zeigte den Weg bzw. den Ausgangspunkt, 
wie ae i Autor im ersten Abschnitt darlegt. 
Die tone tag evidenten Grundtatsachen des menschlichen Verhaltens werden als 
Grundlage der ganzen Theorie analysiert, und zwar wesentlich-tiefer schiirfend 
on von den Begriindern und Nachfahren der Wiener und der Leamater © Schule. 


wird der Ablauf der Wirtschaft unter Verwendung einer még- 
Rebs einfachen Mathematik untersucht, wobei die unvermeidlichen  idealisierenden 
Vora in vieler Hinsicht weniger einschrankend sind als etwa bei der 
Lausanner und ihren Nachfahren, Fiir eine Weiterentwicklung der Theorie 


unter schrittweiser aes dieser idealisierenden Voraussetzungen bildet eine 
Arbeit wie die vorliegende die unentbehrliche Grundlage, da nur so die begriff- 
lichen Zusammenhdnge erstmals geklirt und angemessene ‘wirtschaftliche Begriffe 
exakt definiert werden kénnen 


Zu Kadlekon durch Ihre Buchhandlung 
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